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Abstract. The formation of large-scale vortices is an intriguing phenomenon in two- 
dimensional turbulence. Such organization is observed in large-scale oceanic or atmo- 
spheric flows, and can be reproduced in laboratory experiments and numerical simula- 
tions. A general explanation of this organization was first proposed by Onsager (1949) 
by considering the statistical mechanics for a set of point vortices in two-dimensional 
hydrodynamics. Similarly, the structure and the organization of stellar systems (glob- 
ular clusters, elliptical galaxies,...) in astrophysics can be understood by developing a 
statistical mechanics for a system of particles in gravitational interaction as initiated 
by Chandrasekhar (1942). These statistical mechanics turn out to be relatively similar 
and present the same difficulties due to the unshielded long-range nature of the in- 
teraction. This analogy concerns not only the equilibrium states, i.e. the formation of 
large-scale structures, but also the relaxation towards equilibrium and the statistics of 
fluctuations. We will discuss these analogies in detail and also point out the specificities 
of each system. 



1 Introduction 

Two-dimensional flows with high Reynolds numbers have the striking property of 
organizing spontaneously into coherent structures (the vortices) which dominate 
the dynamics |Q (see Fig. |^). The robustness of Jupiter's Great Red Spot, a 
huge vortex persisting for more than three centuries in a turbulent shear between 
two zonal jets, is probably related to this general property. Some other coherent 
structures like dipoles (pairs of cyclone/anticyclone) and sometimes tripoles have 
been found in atmospheric or oceanic systems and can persist during several 
days or weeks responsible for atmospheric blocking. Some astrophysicists invoke 
the existence of organized vortices in the gaseous component of disk galaxies in 
relation with the emission of spiral density waves . It has also been proposed 
that planetary formation might have begun inside persistent gaseous vortices 
born out of the protoplanetary nebula 121 15,3C,B3 (see Fig. H). As a result, 



hydrodynamical vortices occur in a wide variety of geophysical or astrophysical 
situations and their robustness demands a general understanding. 

Similarly, it is striking to observe that self-gravitating systems follow a kind 
of organization despite the diversity of their initial conditions and their en- 
vironement (see Fig. ^). This organization is illustrated by morphological 
classification schemes such as the Hubble sequence for galaxies and by simple 
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rules which govern the structure of individual self-gravitating systems. For ex- 
ample, elliptical galaxies display a quasi-universal luminosity profile described 
by de Vaucouleur's i?^/^ law and most of globular clusters are well-fitted by 
the Michie-King model. On the other hand, the flat rotation curves of spiral 
galaxies can be explained by the presence of a dark matter halo with a density 
profile decreasing as at large distances. The fractal nature of the interstellar 
medium and the large scale structures of the universe also display some form of 
organization. 




The question that naturally emerges is what determines the particular con- 
figuration to which a self-gravitating system or a large-scale vortex settles. It 
is possible that their actual configuration crucially depends on the conditions 
that prevail at their birth and on the details of their evolution. However, in view 
of their apparent regularity, it is tempting to investigate whether their organi- 
zation can be favoured by some fundamental physical principles like those of 
thermodynamics and statistical physics. We ask therefore if the actual states of 
self-gravitating systems in the universe and coherent vortices in two-dimensional 
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Fig. 2. A scenario of planet formation inside large-scale vortices presumably present 
in the Keplerian gaseous disk surrounding a star at its birth. Starting from a random 
vorticity field, a series of anticyclonic vortices appears spontaneously (upper panel). 
Due to the Coriolis force and to the friction with the gas, these vortices can efficiently 
trap dust particles passing nearby (lower pannel). The local increase of dust concen- 
tration inside the vortices can initiate the formation of planetesimals and planets by 
gravitational instability. This numerical simulation is taken from jl^ ]. 



turbulent flows are not simply more probable than any other possible configura- 
tion, i.e. if they cannot be considered as maximum entropy states. This statistical 
mechanics approach has been initiated by Onsager |101| for a system of point 
vortices and by Chandrasekhar in the case of self-gravitating systems. 

It turns out that the statistical mechanics of two-dimensional vortices and 
self-gravitating systems present a deep analogy despite the very different physical 
nature of these systems. This analogy was pointed out by Chavanis in |29,3^,^ 
and further developed in ||,|^,||j36[|^,|§ . In the following, we will essentially 
discuss the statistical mechanics of 2D vortices and refer to the review of Pad- 
manabhan [[103| (and his contribution in this book) for more details about the 
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Fig. 3. Large-scale structures in the universe as observed with the Hubble space tele- 
scope. The analogy with Fig. hi is striking and will be discussed in detail in this paper. 



statistical mechanics of self-gravitating systems. We will see that the analogy be- 
tween two-dimensional vortices and (three-dimensional) self-gravitating systems 
concerns not only the prediction of the equilibrium state, i.e. the formation of 
large-scale structures, but also the statistics of fluctuations and the relaxation 
towards equilibrium. 

This paper is organized as follows. In Sec. |^, we discuss the statistical me- 
chanics of point vortices introduced by Onsager |101| and further developed by 
Joyce & Montgomery and Pointin & Lundgren |107[ among others (see a 
complete list of references in the book of Newton ||9^ ) . We discuss the existence 
of a thermodynamic limit in Sec. 2/7 and make the connexion with field theory. 
Statistical equilibrium states of axisymmetric flows are obtained analytically in 



Sec. p.8[ - |2.9| . The r elation with equilibrium states of self-gravitating systems is 
shown in Sec. ^.10 In Sec. we discuss the statistics of velocity fluctuations pro- 
duced by a random distribution of point vortices and use this stochastic approach 
to obtain an estimate of the diffusion coefficient of point vortices. Application 
to 2D decaying turbulence is considered in Sec. 3^. In Sec. ^, we describe the 
relaxation of a point vortex in a thermal bath and analyze this relaxation in 
terms of a Fokker-Planck equation involving a diffusion and a drift. In Sec. ^, we 
develop a more general kinetic theory of point vortices. A new kinetic equation 
is obtained which satisfies all conservation laws of the point vortex system and 
increases the Boltzmann entropy (H-theorem). We mention the connexion with 
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the kinetic theory of stars developed by Chandrasekhar In Sec. |[ we dis- 
cuss the violent relaxation of 2D vortices and stellar systems. We mention the 
analogy between the Vlasov and the Euler equations and between the statistical 
approach developed by Lynden-Bell [M for coUisionless stellar systems and by 



Kuz'min p3|, Miller [p5| and Robert & Sommeria |111| for continuous vorticity 



fields. The concepts of "chaotic mixing" and "incomplete relaxation" are dis- 



cussed in the light of a relaxation theory in Sec. 3.3, Application of statistical 



mec hani cs to geophysical flows and Jupiter's Great Red Spot are evocated in 



Sec. 6.4. 



2 Statistical mechanics of point vortices in 
two-dimensional hydrodynamics 

2.1 Two-dimensional perfect flows 

The equations governing the dynamics of an invisicid flow are the equation of 
continuity and the Euler equation: 

^+V{pn) = 0, (1) 

du , ^ 1 

- + (uV)u^--Vp. (2) 

For an incompressible flow, the equation of continuity reduces to the condition 

Vu = 0. (3) 

If, in addition, the flow is two-dimensional, this last equation can be written 
dxU + dyV = 0, where (u, v) are the components of the velocity. According to the 
Schwarz theorem, there exists a streamfunction ip such that u — dytjj, v — ~dxip, 
or, equivalently 

u = -z X VtA, (4) 
where z is a unit vector normal to the flow. The vorticity 

w = V X u = tjz, with Lo = dxV ~ dyU, (5) 

is directed along the vertical axis. According to the Stokes formula, the circu- 
lation of the velocity along a closed curve (C) delimiting a domain area {S) 
is 

r= (f ud\= [ ujd?Y. (6) 

J(C) J(S) 

Taking the curl of Eq. (|^) , we find that the vorticity is related to the stream 
function by a Poisson equation 



Alp = -LU, 



(7) 
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where A = d^^ + dyy is the Laplacian operator. In an unbounded domain, this 
equation can be written in integral form as 

V'(r,i) = -^ J u;{r',t)ln\r-r'\d\', (8) 

and the velocity field can be expressed in terms of the vorticity as 

u(r,t) = ^zx I u^ir', t)^^d\'. (9) 

In a bounded domain, Eq. (||) must be modified so as to take into account 
vortex images. The impermeability condition implies that -0 is constant on the 
boundary and we shall take ■)/' = by convention. Taking the curl of Eq. (|^), the 
pressure term disapears and the Euler equation becomes 

^ + uVcj = 0. (10) 
ot 

This corresponds to the transport of the vorticity uj by the velocity field u. It is 
easy to show that the flow conserves the kinetic energy 

E-J y rf^r. (11) 

Using Eqs. (^(0), one has successively 

E^IJ {^^fd^v - i y ^(-Z\V)d2r = 1 y ^V-rf'r, (12) 

where the second equality is obtained by a part integration with the condition 
'i/' = on the boundary. Therefore, E can be interpreted either as the kinetic 
energy of the flow (see Eq. (pi])) or as a potential energy of interaction between 
vortices (see Eq. (Ilj 



2.2 The point vortex gas 

We shall consider the situation in which the velocity is created by a collection 
of N point vortices. In that case, the vorticity field can be expressed as a sum 
of (5-functions in the form 

N 

a;(r,i)=^7,<5(r-r,(i)), (13) 

where {t) denotes the position of point vortex i at time t and 7^ is its circulation. 
According to Eqs. (p[) (p^), the velocity of a point vortex is equal to the sum of 
the velocities V(j i) produced by the — 1 other vortices, i.e. 

V. = ^itl^ V(j^») = -j^zx i'"^"'''' . (14) 
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As emphasized by KirchhofF , the above dynamics can be cast in a Hamil- 
tonian form 

dt ~ % ' dt ~ dx, ' ^ ^ 

where the coordinates {x,y) of the point vortices are canonically conjugate. 
These equations of motion still apply when the fluid is restrained by boundaries, 
in which case the Hamiltonian ( p^ is modified so as to allow for vortex images, 
and may be constructed in terms of Green's functions depending on the geometry 
of the domain. Since H is not explicitly time dependant, it is a constant of 
the motion and it represents the "potential" energy of the point vortices. The 
other conserved quantities are the angular momentum ^ivf and the impulse 
"Y^^^iVi. Note that the Hamiltonian (|l6| ) does not involve a "kinetic" energy 

2 

of the point vortices in the usual sense (i.e., a quadratic term ^). This 
is related to the particular circumstance that a point vortex is not a material 
particle. Indeed, an isolated vortex remains at rest contrary to a material particle 
which has a rectilinear motion due to its inertia. Point vortices form therefore 
a very peculiar Hamiltonian system. Note also that the Hamiltonian of point 
vortices can be either positive or negative (in the case of vortices of different 
signs) whereas the kinetic energy of the flow is necessarily positive. This is clearly 
a drawback of the point vortex model. 



2.3 The microcanonical approach of Onsager (1949) 



The statistical mechanics of point vortices was first considered by Onsager |101| 
who showed the existence of negative temperature states at which point vortices 
of the same sign cluster into "supervortices" . He could therefore explain the 
formation of large, isolated vortices in nature. This was a remarkable anticipation 
since observations were very scarce at that time. 

Let us consider a liquid enclosed by a boundary, so that the vortices are 
confined to an area A. Since the coordinates {x, y) of the point vortices are 
canonically conjugate, the phase space coincides with the configuration space 
and is finite: 

J dxidyi.-.dxiqdyN = (^j dxdy^ — A^ . (17) 

This striking property contrasts with most classical Hamiltonian systems con- 
sidered in statistical mechanics which have unbounded phase spaces due to the 
presence of a kinetic term in the Hamiltonian. 

As is usual in the microcanonical description of a system of N particles, we 
introduce the density of states 



g{E) = J dxidyi...dxNdyNs(^E ~ H{xi,yi, . 



XN,yN)], (18) 
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which gives the phase space volume per unit interaction energy E. The equiUb- 

rium A^-body distribution of the system, satisfying the normahzation condition 
J IJ,{ri,...,rN)(fri...(frN = 1, is given by 

M(ri, ...,r;v) = -^J{E - H{ru ...,r;v)). (19) 

This formula simply means that, in the microcanonical ensemble, all accessi- 
ble microstates (having the required energy E) are equiprobable at statistical 
equilibrium. 

The phase space volume which corresponds to energies if(ri, rjv) less than 
a given value E can be written 



${E)= / g{E)dE. 



(20) 



It increases monotonically from zero to when E goes from Emin to +oo. 
Therefore, g{E) = d^{E)/dE will have a maximum value at some E = Em, say, 
before decreasing to zero when E +oo. 

In the microcanonical ensemble, the entropy and the temperature are defined 

by 

S = \ng{E), /3 = ^ = ^- (21) 

For E > E„i, S{E) is a decreasing function of energy and consequently the 
temperature is negative. Now, high energy states E ^ E^ are clearly those in 
which the vortices of the same sign are crowded as close together as possible. For 
energies only slighlty greater than Em, the concentration will not be so dramatic 
but there will be a tendency for the vortices to group themselves together on 
a macroscopic scale and form "clusters" or "supervortices" . By contrast, for 
E < E„i, the temperature is positive and the vortices have the tendency to 
accumulate on the boundary of the domain in order to decrease their energy. For 
a system with positive and negative vortices, the negative temperature states, 
achieved for relatively high energies, consist of two large counter- rotating vortices 
physically well separated in the box. On the contrary when E — oo, the 
temperature is positive and vortices of opposite circulation tend to pair off. 



2.4 The equation of state 

For a two-dimensional gas of particles interacting via a Coulombian potential in 
Inr, the equation of state can be derived exactly. Let us assume that the system 
is enclosed in a domain of surface V = E? . The density of states can be written 

g{E,V)=j^ ...J^ lld'r,S^E+-Y,iajl^h-rj\j. (22) 
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Making the change of variable x = r/i?, we find that g{E,V) = g{E' ,1) 
with E' = E + I'aV J2i<j lilj ■ Therefore, the entropy satisfies S{E,V) = 
NlnV + S{E', 1) and the pressure P = T{dS/dV)E is exactly given by 



^^^(l+4i^E^.O,), (23) 

^ i<'J ^ 



If the vortices have the same circulation 7, we obtain 



In particular the pressure vanishes for 



We shall see in Sec. 2^ that this negative critical inverse temperature is the 
minimum inverse temperature that the system can achieve. 

If, on the other hand, we consider a neutral system consisting of N/2 vortices 
of circulation 7 and N/2 vortices of circulation —7, we find 

p=Afi_/!2!y (26) 



This result is well-known is plasma physics |113|. The critical temperature at 
which the pressure vanishes is now positive 

Pc^%, (27) 
7 

and independant on the number of point vortices in the system. For [3 > [3c the 
pressure is negative so this range of temperatures is forbidden. For simplicity, we 
have not taken into account the contribution of images in the previous calcula- 
tions; this can slightly change the results but this should not alter the existence 
of the critical inverse temperatures reported above. 



2.5 The mean-field approximation 

It is easy to show that the exact distribution of point vortices ( [T^ ) expressed 
in terms of (5- functions is solution of the Euler equation (p^. This is proved as 
follows. Taking the derivative of Eq. (Oh with respect to time, we obtain 



— = -^7V5(r-r,(t))V,. (28) 



dt 

z— i 

Since = u{r i{t),t), we can rewrite the foregoing equation in the form 



dt 



--V^7^(r-r,(t))u(r,i). (29) 



10 Pierre-Henri Chavanis 



Since the velocity is divergenceless, we obtain 

= -u(r, t)V ^S{r - r,{t)) = -uVu;. (30) 

i—l 

Therefore, in the point vortex model, the Euler equation (nOh contains exactly 



the same information as the Hamiltonian system (15). 

This description in terms of (S-functions, while being technically correct, is 
useless for practical purposes because it requires the knowledge of the exact tra- 
jectories of the point vortices for an arbitrary initial condition or the solution of 
the Euler equation (p^). When A'' is large, this task is impossibly difficult. There- 
fore, instead of the exact vorticity field expressed in terms of ^-functions, one is 
more interested by functions which are smooth. For that reason, we introduce a 
smooth vorticity field (w) (r, t) which is proportional to the average number of 
vortices contained in the cell (r, r -I- dr) at time t. This description requires that 
it is possible to divide the domain in a large number of cells in such a way that 
each cell is (a) large enough to contain a macroscopic number of point vortices 
but (b) small enough for all the particles in the cell can be assumed to possess 
the same average characteristics of the cell. 

Formally, the average vorticity field is given by 

N 

(c.)(r,t)=^7(^(r-r,(t))), (31) 
1=1 

where the statistical average of a function X{ri, r^r) is defined by 

(X}= J ^l{rl,...,rN,t)X{rl,...,rN)d^rl...d^rN, (32) 

where /i(ri, ...,rAr,t) is the A'^-body distribution function of the system at time 
t. The average vorticity can be rewritten 

(a;)(r,i) -7V7F(r,t) =7(n)(r,t), (33) 

where we have introduced the one-vortex distribution function 

F(ri,t) = J ^i{ri,...,rN,t)d^r2...d^rN, (34) 

and the local vortex density (n) — NP{r, t). In the foregoing relations, we have 
implicitly used the fact that the vortices are identical. The average energy of the 
system is given by 



£;=(ff)=-^^7'(ln|r,:-r,|) 

= — -N{N - 1)7^ / g{r, r', t) In |r - r'\d^rd^r', (35) 
47r J 
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where 



5(ri,r2,i)= J fi{ri, ...,rN,t)(fr2...(frN, (36) 
is the two-body distribution function. In the mean-field approximation, which is 



exact in a properly defined thermodynamic limit with N +00 (see Sec. 2.7), 
we have 

.9(ri,r2,t) = F(ri,t)P(r2,i). (37) 
Accounting that N{N — 1) ~ N'^ for large N, the average energy takes the form 

E = -— f(uj){r,t){Lu){r',t)ln\r-r'\d^rd^r'. (38) 
An J 

Using Eq. (||), the expression ( ^ for E can be rewritten 

E=IJ (^)^d^r = J (39) 

where ip is the streamfunction created by the average vorticity (lu) and where 
(u) is the smooth velocity field. 



2.6 The maximum entropy state 

We now wish to determine the equilibrium distribution of vortices following a 
statistical mechanics approach. Using Boltzmann procedure, we divide the do- 
main in macrocells with area A. Let rii denote the number of point vortices 
in the cell Ai. We now decompose each macrocell into h' microcells with equal 
area. A macrostate is determined by the distribution {rij}. Several configura- 
tions can lead to the same macrostate: each of them will be called a microstate. 
Using a combinatorial analysis, the number of microstates corresponding to the 
macrostate {ui} is 

W{{n,})^Nll[—. (40) 

I 

The logarithm of this number defines the Boltzmann entropy. Using Stirling 
formula and considering the continuum limit in which A^v — > 0, we get the 
classical formula 



-N 



J P(r)lnP(r)d2r, (41) 



where P{r) is the density probability that a point vortex be found in the sur- 
face element centered on r. At equilibrium, the system is in the most probable 
macroscopic state, i.e. the state that is the most represented at the microscopic 
level. This optimal state is obtained by maximizing the Boltzmann entropy ( ^l| ) 
at fixed energy ( p9| ) and vortex number N, or total circulation 

r = Nj= I (Lj)(fr. (42) 
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Writing the variational principle in the form 

SS - f36E - aSr = 0, (43) 

where /3 and a are Lagrange multipliers, it is readily found that the maximum 
entropy state corresponds to the Boltzmann distribution 

(lu) = Ae-I^-'^, (44) 

with inverse temperature p. We can account for the conservation of angular 
momentum L — J{uj)r'^(Pr (in a circular domain) and impulse P ~ J{Lu)y(Pr 
(in a channel) by introducing appropriate Lagrange multipliers f2 and U for 
each of these constraints. In that case, Eq. (^) remains valid provided that we 
replace the streamfunction ip by the relative streamfunction ip' — ^ + — Uy. 
Substituting the Boltzmann relation between {lo) and V' in the Poisson equation 
(^, we obtain a differential equation for the streamfunction 

- A^j = Ae"''^'^, (45) 

which determines the statistical equilibrium distribution of vortices. This Boltz- 
mann Poisson equation can be easily generalized to include a spectrum of circu- 
lations among the vortices. 

The combinatorial analysis presented in this section was performed by Joyce 
& Montgomery JZOt. The mean- field equation (MS) was also obtained by Pointin 



& Lundgren |107| from the equilibrium BGK hierarchy. In fact, this mean-field 
approach was first developed by Onsager but his results were not published (U. 
Frisch, private communication). 



2.7 Field theory and thermodynamic limit 

Let us consider a collection of N point vortices with equal circulation 7 in a 
bounded domain of size R. We want to give a rational to the mean-field approxi- 
mation considered previously by defining a proper thermodynamic limit for point 
vortices. Rigorous results have been established in |p7| , ^ . In the following, we 
present a le ss rigor ous, albeit equivalent, field theory approach inspired by the 
work of ||6^ , 103| , 57| for self-gravitating systems. 



In the microcanonical ensemble, the quantity of fundamental interest is the 
density of states 



g{E)= j 5(^E-H{v^,...,VM)^X{(fv,, 



(46) 



which is explicitly given by 
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where the potential of interaction has been normalized by R. For simplicity, 
we have ignored the contribution of the images but this shall not affect the 
final results. We now introduce the change of variables x = v/R and define the 
function ^ 

m(xi,...,Xa,) = — ^ln|x,: -Xjl, (48) 

and the dimensionless energy 

27r£: 

In terms of these quantities, the density of states can be rewritten 

The proper thermodynamic limit for a system of point vortices with equal cir- 
culation in the microcanonical ensemble is such that N — > +00 with fixed A. We 
see that the box size R does not enter in the normalized energy A. Therefore, the 
thermodynamic limit corresponds to iV ^ +cx) with 7 ^ A^^^ and E ^ \. 
This is a very unusual thermodynamic limit due to the non-extensivity of the 
system. Note that the total circulation F = remains fixed in this process. 
For sufficiently large A^, the density of states can be written 



with 



1 Ppe^^M 8{a^E\p\)\^6{\- ^ KOd^r), (51) 

S[p]=- I pir)\npir)d\ (52) 

E[p]=-^J p{r)p{r')\n\r-r'\d^rd^r'. (53) 

In the above formula, g{E) has been expressed as a functional integral over the 
macrostates recalling that e^'^'^'l = W^({p}) give s the number of microstates 
corresponding to the macrostate p(r) (see Sec. ^.6[ ). The crucial point to realize 
is that the vortex number N appears explicitly in the exponential, all other 
terms being of order unity. Therefore, at the thermodynamic limit N — > +00, 
the functional integral is dominated by the distribution p*(r) which maximizes 
the Boltzmann entropy ( p2[ ) under the constraints of fixed circulation and energy 
brought by the (5-functions. In this limit, the mean-field approximation is exact 
and we have 

g(i?) = e^«[''*l, SiE)^NS[p4. (54) 
In the canonical ensemble, the object of interest is the partition function 

N 
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which is the normahzation factor in the Gibbs measure 



M(ri,...,r^) = |e-^«('^--'^-). (56) 



The free energy is defined by the relation 



F(/3) = -ilnZ. (57) 



Using the notations introduced previously, we can rewrite the integral (jS^) in 
the form 

Z{(3) = V^ ... TTfi^x.e''"^^^'-'''"), (58) 
Jo Jo 

where 

V = -T, — : (59) 

is the normalized inverse temperature. It plays the role of the "plasma param- 
eter" in plasma physics. It can be shown that the partition function (|5^) is 
convergent only for rj > rjc — —4:N/{N — 1) which is consistent with the re- 
sult ( |25[ ) . The proper thermodynamic limit for a canonical distribution of point 
vortices is such that N — -|-cxi with fixed rj. In this limit, 7 ~ and (3 ^ N . 

The partition function is related to the density of states by the Laplace 
transform 

/+00 
dE g{E)e-P^. (60) 
-00 

Therefore, for large N one has, using Eq. (^T[), 



(61) 



with 



J[p] = S[p] - 7jE[p]. (62) 



At the thermodynamic limit, the functional integral is dominated by the distri- 
bution /9*(r) which maximizes the Massieu function J[p] under the constraint of 
a fixed circulation. In this limit, the mean- field approximation is exact and we 
have 

Z(/3) = e^-^I^*! , F = E-TS. (63) 

It can be noted that the critical points of entropy at fixed energy and circulation 
(microcanonical description) and the critical points of free energy at fixed tem- 
perature and circulation (canonical description) coincide. This is not necessarily 
the case for the second order variations of entropy and free energy. Therefore, a 
distribution p(r) can be a maximum of S[p] but a minimum (or a saddle point) 
of J[p\. In that case, the equilibrium is stable in the microcanonical ensemble 
but unstable in the canonical ensemble (and the procedure leading to Eq. ( |6^ ) is 
clearly not valid). When this situation happens, the ensembles are inequivalent 
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and phase transitions occur. This is the case, in particular, for the gravitational 
problem (see Sec. 2.10). 

Finally, the grand canonical partition function is defined by 



+ 00 N r ^ 

where z is the fugacity. Using the Hubbard-Stratanovich transformation 

where G is the Green function of the Laplacian operator, we can rewrite the 
Boltzmann factor in the form 



J e-/^(^«''^'--+V^S".^«('-'), (66) 



where we have assumed /3 < which is the case of physical interest. Substituting 
this result in Eq. (^) and introducing the notation 0(r) = -\/— /37^(r), we can 
easily carry out the summation on N to obtain 

Zac =jv<P e-^/'^^r{l(V^)=^.= e^r)}^ ^^^^ 

Teff = = -zPl^- (68) 

Therefore, the grand partition function of the point vortex gas corresponds to a 
Liouville field theory with an action 

S[<t>\ -7f^ [ rf'r{i(V0)2 - M'e*(r)}. (69) 

While the previous description is formally correct if we define Z and Zqc by 
Eqs. (^5|) and (|6^), it must be noted however that the canonical and grand canon- 
ical ensembles may not have a physical meaning for point vortices. In particular, 
it is not clear how one can impose a thermal bath at negative temperature. 
On the other hand, the usual procedure to derive the canonical ensemble from 
the microcanonical ensemble rests on a condition of additivity which is clearly 
lacking in the present case. 



2.8 Axisymmetric equilibrium states in a disk 



Let us consider a collection of 'N point vortices with circulation 7 confined within 
a disk of radius i?. At statistical equilibrium, the streamfunction -0 is solution 
of the Boltzmann-Poisson equation (45). If we work in a circular domain, we 
must in principle account for the conservation of angular momentum. This can 
lead to bifurcations between axisymmetric and off-axis solutions [117|. We shall, 
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however, ignore this constraint for the moment in order to obtain analytical ex- 
pressions for the thermodynamical parameters. This is a sufficient approximation 
to illustrate the structure of the problem, which is our main concern here. 

If we confine our attention to axisymmetric solutions, the Boltzmann-Poisson 
equation (jiq) can be written 



1 d / 
r dr \ dr 



woe-'^^f'^-'^''). (70) 



i;'{0) = 0, i>{R) - 0, (71) 

where ipo and luq are the values of streamfunction and vorticity at r = 0. In- 
troducing the function (f> = Pj{ip — V'o) and the dimensionless radial distance 
C = {\l3\juJoy^'^r, Eq. (|70| ) can be reduced to the form 

m = '/''(O) = 0, (73) 

with A = 1 if /3 < and X — —1 ii (3 > 0. It turns out that this equation can 
be solved analytically as noticed by a number of authors. With the change of 
variables i = In^ and = 21n^ — z, Eq. (fz^ ) can be rewritten 

This corresponds to the motion of a ficticious particle in a potential V{z) = Ae^. 
This equation is readily integrated and, returning to original variables, we finally 
obtain 

"'^(TTW- '''' 

From the circulation theorem (^ applied to an axisymmetric flow, we have 

dr znr 

where r{r) — uj{r')27rr'dr' is the circulation within r. Taking r = R and 
introducing the dimensionless variables previously defined, we obtain 

7?^^ = -«(/.'(«), (77) 

where a = {\P\"fU!o)^^'^R is the value of ^ at the box radius. Equation ( [tT] ) relates 
a to the inverse temperature rj. Using Eqs. ( [75| ) and ([77|), the vorticity field can 
be written explicitly 
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At positive temperatures (t^ > 0), the vorticity is an increasing function of the 
distance and the vortices tend to accumulate on the boundary of the domain (Fig. 

On the contrary at negative temperatures (77 < 0), the vorticity is a decreasing 
function of the distance and the vortices tend to group themselves in the core of 
the domain to form a "supervortex" (Fig. ||). These results are consistent with 
Onsager's prediction |101|. We also confirm that statistical equilibrium states 
only exist for 77 > r/c = —4, as previously discussed. At this critical temperature, 
the central vorticity becomes infinite and the solution tends to a Dirac peak: 

{u;){r) ^ rSir), for 77 ^ r/, = -4. (79) 
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Fig. 4. Statistical equilibrium states of point vortices at positive temperatures {uj, 
r/nR^). The vortices are preferentially localized near the wall. 



The energy defined by Eq. ( |39| ) can be written in the dimensionless form 
. 2nE 1 
V 



r2 



(80) 



The integral can be carried out explicitly using Eq. ([T^) . Eliminating a between 
Eqs. (80) and ([77|), we find that the temperature is related to the energy by the 
equation of state 



A = 



■In 



V 



1 



(81) 



The corresponding equilibrium phase diagram is represented in Fig. ^ We check 
explicitly that the energy becomes infinite when r] rj^ = —4. The transition 
between positive and negative temperatures occurs for Aq = 1/8. 
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Fig. 5. Statistical equilibrium states of point vortices at negative temperatures showing 
a clustering. For rj = —4, the vortices collapse at the center of the domain and the 
vorticity profile is a Dirac peak. 




Fig. 6. Equilibrium phase diagram (caloric curve) for point vortices with equal circula- 
tion confined within a disk. For simplicity, the angular momentum has not been taken 
into account (i? = 0). 
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Fig. 7. Entropy vs energy plot for a system of point vortices with equal circulation 
confined within a disk. 

The entropy (^l|) can also be calculated easily from the above results. Within 
an unimportant additive constant, it is given by 

f --ln4- fl + -)ln(4 + r;). (82) 

Using Eq. (|8l|), we can easily check that dS/dE — /3, as it should. For A < Aq = 
1/8, the entropy is an increasing function of energy (positive temperatures) and 
for yl > ylo, it is a decreasing function of energy (negative temperatures). For 

A^Q, S/N ~ Inyl and for A +oo, S/N AA (Fig. |). 

It is amusing to note that if we define a local pressure by the relation p(r) = 
(w)/7/3 (which is similar to the local equation of state p = pT for an ideal gas of 
material particles), we find that the pressure P at the boundary of the domain 
is exactly given by Eq. (p^). Furthermore, using the Boltzmann distribution 
(p4[), we easily check that Vp = —{u!)Vtp, which is similar to the equation of 
hydrostatic equilibrium for a fluid in a gravitational field. It is not clear whether 
these results bear more significance than is apparent at first sight. 

2.9 Equilibrium states in an unbounded domain 

In an infinite domain, it is necessary to take into account the conservation of 
angular momentum, because this constraint determines the typical size of the 
system. Furthermore, we shall assume that the vortices have the same sign oth- 
erwise positive and negative vortices will form dipoles and escape to infinity. 
Therefore, there is no equilibrium state in an infinite domain for a neutral sys- 
tem of point vortices. When the conservation of angular momentum is taken into 
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account, the density of states is given by 

giE,L)= I s(E-H{ri,...,rN)) s(L-J2irf)Y[d^r„ (83) 

and the angular velocity of the flow by 

« = 2r(|f) (84) 



Using the same trick as in Sec. 2.4 with the change of variables x = r/\/L, one 
finds the exact result 

-|^(-^)- (-) 

Therefore, the vorticity field is determined by the Boltzmann distribution 

(lu) = Ae-'^^'f'' , (86) 
where ip' is the relative streamfunction 

^' = ^ + ^r^ ^^+^{A + r]y. (87) 

For 77 = 0, one has 

^''^ = ^' ' • 

For large r, the asymptotic behavior of Eq. ( |8^ ) is 

(w) - ;^e-*(^+'')'-' (r ^ +00), (89) 

where we have used V' ~ ^(r/'^Tr) Inr at large distances. From Eq. (p9|), one sees 
that 77 > —4 is required for the existence of an integrable solution. 
Inserting the relation ( |8^ ) in the Poisson equation (^, we get 

- All)' = Ae"—''' -i.'^ + ri). (90) 

2'KLr] 

With the change of variables 



-iNV''^ ( la \ 27r77 



the Boltzmann-Poisson equation ( |90| ) can be written 



(91) 



+ 7^^ = 2^,76* -(4 + ??), (92) 
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and the vorticity ( |8^ becomes 

(c) = -^e*. (93) 

Equation ( |9^ ) has been obtained by Lundgren & Pointin from the equihb- 
rium BGK hierarchy. In the hmit rj — > +00, we have 

(c.)^^, if r<{2L/rY'\ (94) 

and {lj) = otherwise. This vortex patch is the state of minimum energy at 
fixed circulation and angular momentum. For rj —4, one has approximately 

The first factor is an exact solution of Eq. (p3) with the second term on the right 



hand side neglected (see Sec. 2.8). The second factor is a correction for large r, 
in agreement with the asymptotic result expressed by Eq. (^) . The parameter A 
tends to infinity as ry ^ —4 and is determined from the condition J{u!)(Pr — F 
by the formula 

ttA + ln(7rA) =-C- ln(l + |) , (96) 

where C — 0.577... is the Euler constant. For rj = —4, the vorticity profile is a 
Dirac peak and the energy tends to +00. 



2.10 The gravitational TV-body problem 

It is interesting to compare the previous results with those obtained in the case 
of self-gravitating systems [ 103 1 . Formally, the structure of the iV- vortex problem 
shares some analogies with the gravitational A^-body problem. The force by unit 
of mass experienced by a star is given by 

r, = ^r(j->z), F{j ^i) = Gm^^l^, (97) 

where F(j' — > i) is the force created by star j on star i. The force can be written 
as the gradient F = — of a gravitational potential (p which is related to the 
stellar density 

N 

p(r,i) =^TO(5(r-r,), (98) 

i=l 



by the Poisson equation 

A-P = iirGp. 



(99) 
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Furthermore, the equations of motion (Newton's equations) can be put in the 
Hamiltonian form 

dvi dH dvi dH 

dt dvi ' dt dr-i ' 

^ z=l t<j "^^1 

In the analogy between steUar systems and two-dimensional vortices, the star 
density p plays the role of the vorticity w, the force F the role of the velocity V 
and the gravitational potential the role of the streamfunction tj). The crucial 
point to realize is that, for the two systems, the interaction is a long-range un- 
shielded Coulombian interaction (in 13 = 3 or 13 = 2 dimensions). This makes 
the connexion between point vortices and stellar systems deeper than between 
point vortices and electric charges for example. In particular, point vortices can 
organize into large scale clusters, like stars in galaxies, while the distribution of 
electric charges in a neutral plasma is uniform. There are, on the other hand, 
fundamental differences between stars and vortices. In particular, a star creates 
an acceleration while a vortex creates a velocity. On the other hand, the gravi- 
tational interaction is attractive and directed along the line joining the particles 
while the interaction between vortices is rotational and perpendicular to the line 
joining the vortices. 

Despite these important differences, the statistical mechanics of 2D vortices 
and stellar systems are relatively similar. Like the point vortex gas, the self- 
gravitating gas is described at statistical equilibrium by the Boltzmann distri- 
bution 

(p) = Ae-^*, (101) 

obtained by maximizing the Boltzmann entropy at fixed mass M and energy E . 
Its structure is therefore determined by solving the Boltzmann-Poisson equation 

A$ = AnGAe-^'^, (102) 

where A and (3 > have to be related to AI and E. This statistical mechanics 
approach has been developed principally for globular clusters relaxing towards 
equilibrium via two-body encounters |0] . It is clear that the Boltzmann-Poisson 



equation (102) is similar to the Boltzmann-Poisson equation ( [45| ) for point vor- 
tices at negative temperatures. The density profile determined by these equations 
is a decreasing function of the distance, which corresponds to a situation of clus- 
tering (see Figs. |^ and H). The similarity of the maximum entropy problem for 



stars and vortices, and the Boltzmann-Poisson equations ( 102 ) (|45|), is a first 
manifestation of the formal analogy existing between these two systems. 

However, due to the different dimension of space {D = 3 for stars instead 
of -D = 2 for vortices), the mathematical problems differ in the details. First of 
all, the density profile determined by the Boltzmann-Poisson equation ( |102| ) in 
r^^ at large distances leading to the so-called infinite mass 
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Fig. 8. Density profile of tlie self-gravitating gas at statistical equilibrium. The dasiied 
line corresponds to the singular solution p = l/27rG/3r^. 



problem since M = J^°° p Airr^dr ^ +oo There is no such problem for 

point vortices in two dimensions: the vorticity decreases like r""^, or even more 
rapidly if the conservation of angular momentum is accounted for, and the total 
circulation F = uj iirrdr is finite. The infinite mass problem implies that 
no statistical equilibrium state exists for open star clusters, even in theory. A 
system of particles in gravitational interaction tends to evaporate so that the 
final state is just two stars in Keplerian orbit. This evaporation process has been 
clearly identified in the case of globular clusters which gradually lose stars to 
the benefit of a neighboring galaxy. In fact, the evaporation is so slow that we 
can consider in a first approximation that the system passes by a succession of 
quasiequilibrium states described by a truncated isothermal distribution function 
(Michie-King model) This justifies the statistical mechanics approach in that 
sense. Another way of solving the infinite mass problem is to confine the system 
within a box of radius R. However, even in that case, the notion of equilibrium 
poses problem regarding what now happens at the center of the configuration. 

The equilibrium phase diagram (i?, T) for bounded self-gravitating systems is 
represented in Fig. |^. The caloric curve has a striking spiral behavior parametrized 
by the density contrast TZ — p{0)/p{R) going from 1 (homogeneous system) to 
-|-oo (singular sphere) as we proceed along the spiral. There is no equilibrium 
state below Ec = -0.335GM^/R or = In that case, the system 

is expected to collapse indefinitely. This is called gravothermal catastrophe in 
the microcanonical ensemble (fixed E) and isothermal collapse in the canonical 
ensemble (fixed T) . Dynamical models show that the collapse is self-similar and 
develops a finite time singularity [|l05[|85| , p6| , ^ 84,4^. However, although the 
central density goes to +oo, the shrinking of the core is so rapid that the core 
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Fig. 9. Equilibrium phase diagram for self-gravitating systems confined within a box. 
For sufficiently low energy or temperature, there is no equilibrium state and the system 
undergoes gravitational collapse. 
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Fig. 10. Equilibrium phase diagram for self-gravitating fermions |^]. The degeneracy 
parameter /i plays the role of a small-scale cut-off e ~ 1/jU. For e ^ 0, the classical 
spiral of Fig. O is recovered. 
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mass goes to zero. Therefore, the singularity contains no mass and this process 
alone cannot lead to a black hole. 

Since the T{E) curve has turning points, this implies that the microcanon- 
ical and canonical ensembles are not equivalent and that phase transitions will 
occur 1 103 1 . In the microcanonical ensemble, the series of equilibria becomes un- 
stable after the first turning point of energy (MCE) corresponding to a density 
contrast of 709. At that point, the solutions pass from local entropy maxima to 
saddle points. In the canonical ensemble, the series of equilibria becomes unsta- 
ble after the first turning point of temperature (CE) corresponding to a density 
contrast of 32.1. At that point, the solutions pass from minima of free energy 
{F ^ E — TS) to saddle points. It can be noted that the region of negative 
specific heats between (CE) and (MCE) is stable in the microcanonical en- 
semble but unstable in the canonical ensemble, as expected on general physical 
grounds. The thermodynamical stability of isothermal spheres can be deduced 
from the topology of the (3 — E curve by using the turning point criterion of 
Katz who has extended Poincare's theory of linear series of equilibria. The 
stability problem can also be reduced to the study of an eigenvalue equation 
associated with the second order variations of entropy or free energy as studied 



by Padmanabhan |102| in the microcanonical ensemble and by Chavanis p37| in 
the canonical ensemble. This study has been recently extended to other statisti- 
cal ensembles grand canonical, grand microcanonical, isobaric... The same 
stability limits as Katz are obtained but this method provides in addition the 
form of the density perturbation profiles that trigger the instability at the critical 
points. It also enables one to show a clear equivalence between thermodynamical 
stability in the canonical ensemble and dynamical stability with respect to the 
Navier-Stokes equations (Jeans problem) ||3^,0. These analytical methods can 
be extended to general relativity |38| . It must be stressed, however, that the 
statistical equilibrium states of self-gravitating systems are at most metastable: 
there is no global maximum of entropy or free energy for a classical system of 
point masses in gravitational interaction 

Phase transitions in self-gravitating systems can be studied in detail by in- 
troducing a small-scale cut-off e in order to regularize the potential. This can 
be achieved for example by considering a system of self-gravitating fermions 
(in which case an effective repulsion i s played by the Pauli exclusion principle) 
32lg,|5^, Pj43 | or a hard spheres gas fl |l03| , |l20| . Other forms of regularization are 



possible 1 59 128, 45|. For these systems, there can still be gravitational collapse 



but the core will cease to shrink when it feels the influence of the cut-off. The 
result is the formation of a compact object with a large mass: a "fermion ball" or 
a hard spheres "condensate" . The equilibrium phase diagram of self-gravitating 
fermions is represented in Fig. |l^ and has been discussed at length by Chavanis 
in the light of an analytical model. The introduction of a small-scale cut- 
off has the effect of unwinding the classical spiral of Fig. ^. For a small cut-off 
e ^ 1, the trace of the spiral is still visible and the T{E) curve is multivalued 
(Fig. ^l|). This can lead to a gravitational first order phase transition between a 
gaseous phase with an almost homogeneous density profile (upper branch) and a 
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Fig. 11. Enlargement of the phase diagram near the tricritical point in the micro- 
canonical ensemble. A priori, the phase transition should occur at the energy -Et(/i) at 
which the gaseous phase and the condensed phase have the same entropy. In fact, the 
entropic barrier played by the unstable solution on the wiggling branch is so hard to 
cross that the transition will not occur at Et but rather at, or near, Ec the point of 
gravothermal catastrophe ^,0). 



condensed phase with a core-halo structure (lower branch). At a critical cut-off 
value e = emtPj the two phases merge and the gravitational phase transition 
disapears. This particular point is sometimes called a microcanonical tricritical 
point (MTP). 

For e > tMTP, the curve has the form of Fig. |l2|. The T{E) curve is now 
univalued so that the equilibrium states are always stable in the microcanonical 
ensemble (they are global entropy maxima). In particular, the region of nega- 
tive specific heats (leading to a convex dip in the entropy vs energy plot |Q ) is 
allowed in the microcanonical ensemble. By contrast, the E(T) curve is multi- 
valued and this can lead to a normal first order phase transition in the canonical 
ensemble. The gaseous and condensed phases are thus connected by a Maxwell 
plateau which replaces the region of negative specific heats. At e = cctp, the 
two phases merge, the specific heat becomes infinite and the phase transition 
is second order. This particular point is sometimes called a canonical tricritical 
point (CTP). For e > ecTP, the T{E) curve is monotonic and the specific heat is 
positive. Therefore, as the cut-off parameter e increases, the self-gravitating gas 
consecutively exhibits gravitational first order, normal first order, second-order 



and no phase transition at all |45 . A similar behavior exists for the Blume- 
Emery-Grifhths (BEG) model with infinite range interaction and is probably 
representative of other systems presenting inequivalence of ensembles (small sys- 
tems or systems with long-range interactions). 
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Fig. 12. Same as Fig. |l^ near the tricritical point in the canonical ensemble. The 
region of negative specific heats, allowed in the microcanonical ensemble, is replaced 
by aphase transition in the canonical ensemble. The notion of metastability is discussed 

in mM- 



The equilibrium phase diagram of point vortices in two dimensions is com- 
pletely different (Fig. |6|). Since there is no turning point of energy or tempera- 
ture, we can immediately infer that the equilibria are always stable and that the 
microcanonical and canonical ensembles are equivalent. In the microcanonical 
ensemble, an equilibrium state exists for all values of energy. Therefore, there 
is no "gravothermal catastrophe" in two dimensions The solutions of the 
Boltzmann-Poisson system for arbitrary dimension D and the disappearan ce of 
the spiral as we approach the critical dimension D — 2 have been studied in | 116 | 
(note that the spiral also disapears for D > 10). There exists, on the other hand, 
a critical inverse temperature rjc = —4 in two dimensions below which there is 
no equili briu m state in the canonical ensemble. This can lead to a situation of 
collapse [ 116 but, as indicated previously, it is not clear whether the canonical 
ensemble makes sense for a gas of point vortices. 



3 Statistics of velocity fluctuations arising from a random 
distribution of point vortices 

3.1 The marginal Gaussian distribution 

The aim of equilibrium statistical mechanics is to predict the final configuration 
of a system resulting from a complex evolution. We seek now to develop a kinetic 
theory of point vortices to determine how the system will reach this equilibrium 
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configuration. The first problem to consider is the characterization of the velocity 
fluctuations experienced by a point vortex. These fluctuations are responsible for 
a diffusion process which is one of the driving source of the evolution. In this 
section, we describe a simple stochastic model | |6^ , |97 ,126 55|, p^ , [i7| . [18[|4C| | from 



which we can obtain an estimate of the diffusion coefficient of point vortices 
[|7| . We follow the presentation of Chavanis & Sire . 

Let us consider a collection of N point vortices randomly distributed in a disk 
of radius R. We assume that the vortices have a Poisson distribution, i.e. their 
positions are independent and uniformly distributed in average over the entire 
domain. In terms of statistical mechanics, this corresponds to an equilibrium 
state with /3 = in a domain with no specific symmetry or /3 — > +00 (i.e. a state 
of minimum energy) if the angular momentum is conserved. We are particularly 
interested in the "thermodynamic limit" in which the number of vortices and 
the size of the domain go to infinity {N 00, i? — > 00) in such a way that the 
vortex density n = remains finite. In this limit, the Poisson distribution is 
stationary and is well-suited to the analysis of the fluctuations. For simplicity, 
we assume that the vortices have the same circulation 7. There is therefore a 
solid rotation of the system (V) = ■^njr±. We shall work in a rotating frame of 
reference so as to ignore this solid rotation. The generalization of our results to 



a spectrum of circulations among the vortices is relatively straightforward 1 47 

The velocity V occurring at the center O of the domain is the sum of the 
velocities $i (i — 1, TV) produced by the vortices: 



N 



27r r. 



2 , (103) 



where denotes the position of the i*'' vortex relative to the point under consid- 
eration and, by definition, r±i is the vector rotated by +t:/2. Since the vortices 
are randomly distributed, the velocity V fluctuates. It is therefore of interest to 
study the statistics of these fluctuations, i.e. the probability W(V)(PY that V 
occurs between V and V + dV. Basically, we have to determine the distribution 
of a sum of random variables. If the variance of the individual velocity $ were 
finite, we could immediately apply the central limit theorem and deduce that the 
distribution of V is Gaussian. Alternatively, if the variance of $ were rapidly 
diverging (i.e., algebraically), the distribution of V would be a Levy law with 
an infinite variance and an algebraic tail. In the present case, the problem is 
intermediate between these two situations because the variance of the velocity 
created by a single vortex 

(^^) . £^ rir)^^d^r . ^^2.rdr, (104) 



diverges slowly, i.e. logarithmically. For that reason, the distribution of V will 
be intermediate between Gaussian and Levy laws. 
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The distribution Wn (V) of the velocity can be expressed as 

N / ^ \ 

Wn(V)= l[Tir,)d^r,siY-J2^r), (105) 
i=i ^ i=i ' 

where T(vi)(Pvi governs the probability of occurrence of the i-th point vortex 
at position r^. In writing this expression, we have assumed that the vortices are 
identical and uncorrelated. Now, using a method originally due to Markov, we 
express the (5-function appearing in Eq. ( |105[ ) in terms of its Fourier transform 



(27r)2 

With this transformation, Wn^) becomes 



'^(^) = IW^ I ^'"^d'p- (106) 



^N(y) = ^J AM{p)e-'P^<fp, (107) 



with 

N 



1=0 

where we have written 

7 r 



(108) 



* = --^^. (109) 
If we now suppose that the vortices are uniformly distributed on average, then 

r(r) = (110) 



and Eq. (108) reduces to 



Since 



^Nip) I e^P^dPv . (Ill) 

r|=0 / 



dPr^l, (112) 



'|r|=0 

we can rewrite our expression for A]\[{p) in the form 

AT 

A^iP)=\^-^ I a-e^)d'r\ . (113) 



We now consider the limit in which the number of vortices and the size of the 
domain go to infinity in such a way that the density remains finite: 

N 

N — !■ oo, R — !■ cx), n = finite. 
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If the integral occurring in Eq. (|113|) increases less rapidly than N, then 

= e-"<^(^\ (114) 



with 



Cifi) = / {l-e:l^)(fY. (115) 

J|r|=0 



We have dropped the subscript N to indicate that the limit ^ oo, in the 
previous sense, has been taken. Note that A{p) can still depend on N through 



logarithmic factors, so that Eq. (114) must be considered as an equivalent of 
Eq. ( 108 ) for large N , not a true limit. In fact, it is appropriate to consider that 
N — > -l-cxD but not In N since in physical situations the typical number of vortices 
does not exceed 10'*. 

The characteristic function C(p) can be calculated explicitly by taking $ 
as a variable of integration instead of r and transforming to polar coordinates 

. The final result remains complicated but the following expression provides 
a sufficient approximation for our purposes: 



Ci,;,ffJ±l^^Y (U6) 



Since C (p) diverges weakly with N (logarithmically) , the limiting process leading 
to formula ( |114| ) is permissible. The velocity distri bution W(V) is simply the 
Fourier transform of A{p) with the expression ( |ll(i| ) for C{p). This leads to the 
following distribution for the velocity fluctuations: 

W{V)^ e-^^""' {V <l^V,„t{N)), (117) 

in iV 

2 

(y^ycr^t{N)), (118) 
/ 2 \ 1/2 

V^r^t{N)^ [^XnN] Xn'/^ilnN). (119) 



For small fluctuations, the core of the distribution is Gaussian as if the central 
limit theorem were applicable (this is due to the quadratic behavior of C{p)). 
For sufficiently large values of V , the velocity distribution VF(V) decays alge- 
braically as for a Levy law (this is due to the logarithmic term in C{p)). In the 
mathematical limit IniV — > -|-oo, the algebraic tail is rejected to infinity and the 
distribution is purely Gaussian. This is consitent with the generalized form of 
the central limit theorem described by Ibragimov & Linnik and used by Min 



et al. |9^ and Weiss et al. [126|. However, the convergence is so slow with N 
that the algebraic tail is always visible for point vortex systems |^,^ 12£ ■ 



This algebraic tail arises because we are on the frontier between Gaussian and 
Levy laws (see Fig. 1.1 of Bouchaud & Georges ||ll|). Therefore, we proposed to 



call Eqs. (117)-(|119|) the Marginal Gaussian distribution 
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Since the distribution V7(V) decreases like V^^ for V ^ oo, the variance 
of the velocity diverges logarithmically. Since V — 7/27rr, this corresponds to 
a divergence at small scales (r ^ 0). On the other hand, if we were to extend 



the Gaussian distribution (117) for all values of V, we would conclude that its 
variance 

2 

{V^) = !^\nN, (120) 

diverges logarithmically when N ^ oo. This corresponds to a divergence at large 
scales (i? +oo). We can recover these results more directly by calculating the 



variance of V from Eq. (103). Indeed, 



{V^)=Nr ^^2^rdr=^lnf-^Y (121) 

where Lmm is a lower cut-off and R an upper cut-off played here by the size of 
the domain. 

We can give a physical interpretation of the algebraic tail of the velocity 
distribution in terms of the nearest neighbor approximation. Let us calculate 
the velocity distribution due to the nearest neighbor. For that purpose, we must 
first determine the probability Tn.n{'f)dr that the position of the nearest neighbor 
occurs between r and r + dr. Clearly, Tn.n{'r)dr is equal to the probability that 
no vortices exist interior to r times the probability that a vortex (any) exists in 
the annulus between r and r + dr. Therefore, it must satisfy an equation of the 
form 

T„.„(r)dr ={^- T,,,n{r')dr'^n2Tirdr, (122) 
where n = denotes the mean density of vortices in the disk. Differentiating 



Eq. (122) with respect to r we obtain 

.n{r) 



d_ 

dr 



2T:nr 



.W- (123) 



This equation is readily integrated with the condition t„ „(r) ~ 27rnr as r — > 0, 
and we find 

T„.„(r) = 2^nre-™'-'. (124) 

This is the distribution of the nearest neighbor in a random distribution of 
particles. From this formula, we can obtain the exact value for the "average 
distance" d between vortices. By definition, 

Tn.n{r)rdr. (125) 



Hence, 



d=-^. (126) 
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If we assume that the velocity V is entirely due to the nearest neighbor, then 

W„.„(V)d2v = T„.„(r)d2r, (127) 

with 

V = (128) 



Using Eq. (124), we obtain 



The nearest neighbor approximation is expected to give relevant results only for 
large values of the velocity. Thus, we can make the additional approximation 

- 4^^, (130) 



in perfect agreement with Eq. (118) valid for V ^ Vcrit- This shows that the 
algebraic tail of the velocity distribution is essentially produced by the nearest 
neighbor as for a Levy law. 

Note also that the typical velocity due to the nearest neighbor is Vn.n ~ 
^/2iTd where d is the inter- vortex distance (|126|). Hence ~ •ri'j^ j-K^ . Com- 



paring this result with Eq. (120), we see that the velocity due to all the vortices 
is, up to a logarithmic factor, of the same order as the velocity due to the nearest 
neighbor. We can say, in some sense, that the velocity fluctuation is marginally 
dominated by the nearest neighbor and that collective effects are responsible for 
logarithmic corrections. In fact, we can show |Q that the "effective velocity" 
created by a vortex at a distance r from the point under consideration is given 

by 

where yl ~ d is of the order of the inter- vortex distance. For r ^ d, the "effective" 
velocity decays as r^"^ instead of the ordinary law recovered for r d. This 
result indicates that, in a statistical sense, the velocity produced by a vortex is 
shielded by "cooperative" effects. 



3.2 The speed of fluctuations 

The function W^(V) does not provide us with all the necessary information con- 
cerning the fluctuations of V. An important aspect of the problem concerns the 
speed of fluctuations, i.e. the typical duration T(y) of the velocity fluctuation 
V. This requires the knowledge of the bivariate probability W{'V, A)d^Yd^A to 
measure simultaneously a velocity V with a rate of change 

dV " 

A = = (132) 
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7 _ 2(r,,v,)r_Lt 

2 ^4 



27r 



(133) 



where = dvi/dt is the velocity of vortex i. Then, the duration T{V) can be 
estimated by the formula 

nv) = J2=, ,134) 

where 

^ " my) ' ^^^^> 

is the mean square acceleration associated with a velocity fluctuation V. 

To determine the speed of fluctuations, we need therefore to calculate the 
bivariate probability WnCV, A) to measure simultaneously a velocity V with a 



rate of change A = dV/dt. According to Eqs. ( |103| ), ( |l32D and ( |l33D , V and 
A are the sum of N random variables $i and depending on the positions 
Fj and on the velocities of the point vortices. Unlike material particles, the 
variables {r^, v^}, for different i's, are not independent because the velocities of 
the vortices are determined by the configuration {vi} of the system as a whole 
(see Eq. ( |103| )). To be able to solve the problem analytically, we shall make 
a decorrelation approximation and treat {r^, v^} (i — 1,...,N) as independent 
variables governed by the distribution 

1 



tt/c^ 717^ in A* 



resulting from Eqs. (110) and (117). It is remarkable that the distribution ( |136 
is formally equivalent to the Maxwell-Boltzmann statistics of material particles 
at equilibrium. Owing to this analogy, we can interpret the variance 

2 

^=I^lniV, (137) 
47r 

as a kind of kinetic "temperature" of the point vortices. It should not be confused 
with the temperature (3 introduced in Sec. |^ and measuring the clustering of the 
vortices (in the present case, /? = since the vortices are uniformly distributed 
in average). 

When this decorrelation hypothesis is implemented, a straightforward gener- 



alization of the method used in Sec. 3.1 yields 



Wn{V,A) = y^J AMp,<T)e-'^P^^^^d^pd'a, (138) 



with 



AN(p,a) = ( r T(r,v)e*(P*+<^)d2j.^2^ 

V^|r|=0 ^|v|=0 



N 



(139) 



where we have defined 
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1 (^^ _ 2(rv)rj 

2 „4 



27r 



(141) 



Using an integration by parts, the conditional moment of the acceleration for a 
given velocity fluctuation V can be expressed as 

1 r f) A 

W{Y){A^)^ ^(p,0)e-^^dV (142) 

We therefore need to Taylor expand the function A{p,a) for \a\ and carry 



out the integration in Eq. (142). The calculation is relatively tricky but can be 
done analytically B^. Substituting the resulting expression for {A^)y in Eq. 



(134), we obtain the following estimate for the speed of fluctuations: 



^(^)- ,3/2;2 1,^ ^"^ iV«V.ru{N)), (143) 

T{V) = -^==^ {V:^V,ru{N)). (144) 

For weak and large fluctuations, T{V) decreases to zero like V and re- 
spectively. These asymptotic behaviors have a clear physical meaning in the 
nearest neighbor approximation. When r = ^/2t:V is small, corresponding to 
large velocities, it is highly improbable that another vortex will enter the disk 
of radius r before long. By contrast, on a short time scale T ^ r /v ^ {'y /v)V~^ , 
the vortex will have left the disk. When r = ^ /2-kV is large, corresponding to 
small velocities, the probability that the vortex will remain alone in the disk 
is low. The characteristic time before another vortex enters the disk varies like 
the inverse of the number of vortices expected to be present in the disk, i.e. 
T ~ {r /v)l / mir'^ ^ {\/n'fv)V . The demarcation between weak and strong fluc- 
tuations corresponds to V ^ 771-'^/^, i.e. to the velocity produced by a vortex 
distant n~^/^ from the point under consideration. These asymptotic behaviors 
are also consistent with the theory of Smoluchowski concerning the persistence 
of fluctuations (see discussion in Q). 

The average duration of a velocity fluctuation is defined by 



(T) = / T{V)W{V)2'KVdV. (145) 

To leading order in In A^, we obtain 

w \ 1/2 1 
4 /tt \ ' 1 



This formula shows that the typical duration of a velocity fiuctuation scales like 

n7\/In~/V ^ ^ 

This corresponds to the typical time needed by a vortex moving with a typical 
velocity (T/2)i/2 ( 

given by Eq. ( |l2Cl| )) to cross t he in terparticle distance d ~ 
n~^/^, as expected from general physical grounds |126|. 
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3.3 The diffusion coefficient 

According to the previous discussion, we can characterize the fluctuations of 
the velocity of a point vortex by two functions: a function M^(V) which governs 
the occurrence of the velocity V and a function T{V) which determines the 
typical time during which the vortex moves with this velocity. Since the velocity 
fluctuates on a typical time Ttyp = d/ yvV^ which is much smaller than the 
dynamical time to = R/ needed by the vortex to cross the entire domain, 

the motion of the vortex will be essentially stochastic. If we denote by P(r, t) 
the probability density that the particle be found in r at time t, then P(r, t) will 
satisfy the diffusion equation 

dP 

^ = DAP. (148) 



If the particle is at r = Tq at time i = 0, the solution of Eq. (148) is clearly 



^(r,<|ro) = ^e-^, (149) 

where D is the diffusion coefflcient. The mean square displacement that the 
particle is expected to suffer during an interval of time At large with respect to 
the fluctuation time Ttyp, is 

{{Arf) = AD At. (150) 

We can obtain another expression for ((Z\r)^) in terms of the functions W{y) 
and T{V) deflned in the previous sections. Indeed, dividing the interval 

/t+At 
V{t')dt', (151) 

into a succession of discrete increments in position with amount T{Vi)'Vi, we 
readily establish that 

{(Arf) = {T{V)V^)At. (152) 



Combining Eqs. (15C) and ([152) we obtain an alternative expression for the 



diffusion coefficient in the form 



D = \j T{V)W{Y)V'^d'^Y. (153) 



Substituting for T{V) and W{y) in the foregoing expression, we obtain to lead- 
ing order in In TV ||47]| , 

1 /6^^/' 



We should not give too much credit to the numerical factor appearing in Eq. 



(154) since the definition (134) of T{y) is just an order of magnitude. Note that 



the scaling form of D is consistent with the expression 



D ~ Ttyp{V^) ^ 7%/hriV, (155) 
that one would expect on general physical grounds. 
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3.4 Application to 2D decaying turbulence 

The previous results have some direct applications to the context of 2D decaying 
turbulence. The relaxation of 2D decaying turbulence is a three-stage process: 
during an initial transient period, the fluid organizes itself from random fluctu- 
ations and a population of coherent vortices emerges. Then, when two like-sign 
vortices come into contact they merge and form a bigger structure. As time goes 
on, the vortex number decreases and their average size increases, in a process 
reminiscent of a coarsening stage. Finally, when only one dipole is left, it de- 
cays diffusively due to inherent viscosity. Direct numerical simulations and 
experiments |6l| show that the typical core vorticity uj remains constant during 
the course of the evolution and that the density decreases like n ~ with 
^ w 0.7. As the energy E ^ J u.^(Pr ~ nw^a^ ^ na'* is conserved throughout the 
evolution, the typical vortex radius increases like a ~ t^/'^. 

In the punctuated Hamiltonian model ]125| , the flow is approximated by a col- 
lection of vortices with constant vorticity ±uj (and circulation 7 ~ ±i^a|) whose 
centers follow the Hamiltonian dynamics (|l^)(|l^). When two like-sign vortices 
with radii ai and 02 enter in collision, the resulting vortex keeps the same vor- 
ticity and takes a radius a such that a** = af -I- a| which ensures the conservation 
of energy (as discussed above) . This effective model reproduces the results of the 
numerical simulations and experiments, with again ^ ~ 0.7. Since the average 
distance between vortices, of order d ^ t^/'^, increases more rapidly than their 
size, the point vortex model should provide increasing accuracy. We emphasize, 
however, that in the above mentioned studies the density only changes by a fac- 
tor of order 4^5 between the initial time and the final time so that the scaling 
exponent is measured on less than one decade. 

The punctuated Hamiltonian model has been re- investigated recently by Sire 



& Chavanis [115| using a renormalization group procedure which allows for much 
longer time simulations that could otherwise be achieved. It is found that the 
scaling regime is achieved for very late times and is characterized by a decay 
exponent ^ = 1 (an effective exponent ^ ~ 0.7 is recovered for shorter times). 
In addition, the decay of the total area occupied by the vortices results in a 
physical process by which merging occurs principally via 3-vortex collisions in- 
volving vortices of different sign. A simple kinetic theory based on an effective 
3-vortex interaction returns this value f = 1. These theoretical results tend to 
be confirmed by recent direct numerical simulations | p6[ . 

During the decay, the vortices diffuse with a coefficient _D ~ 7 given by Eq. 



(154), where 7 ^ uoa is their circulation (we ignore here logarithmic corrections). 



If the diffusion coefficient were constant, then the dispersion of the vortices 

(r^) - Dt, (156) 

would increase linearly with time as in ordinary Brownian motion. However, 
since D varies with time according to 

D ^ujo^ (157) 



2D vortices and stellar systems 



37 



we expect anomalous diffusion, i.e. 



(158) 



with V ^ \. Substituting Eq. (157) in Eq. (156), we obtain the following relation 
between v and ^: 

^. = 1 + 1 (159) 
This formula is in perfect agreement with long time numerical simulations [v = 



3/2 for f = 1) [115| and experiments {vexp ~ 1-3 — 1.4 in the regime where 
^ = 0.7)pl|. This hyperdiffusive behavior can be interpreted in terms of Levy 



flights |115|, with a large time flight distribution P(t) ~ t ^ with /i = 3 — ^/2 



in agreement with experiments [iiexp ^ 2.6 ± 0.2 for ^ = 0.7) 1 61 



3.5 The spatial correlations in the velocities arising from a random 
distribution of point vortices 



The stochastic approach described in Sec. 3.1 can be generalized to obtain exact 
results concerning the spatial correlations of the velocity fluctuations |Q . First 
of all, the distribution of the velocity increments between two neighboring points 
is given by the 2D Cauchy law §|,|3,|8): 



4|^Vp 



which is a particular Levy law. It is also possible to determine an analytical 
expression for the conditional moment ((5V)v- We find: 

(<SV)v=-^l''|irx+2^^^v| (V»V„„{N)), (162) 
where B{x) denotes the function 

B(x) = i(e^-l-x). (163) 

X 

On the other hand, the spatial auto-correlation function of the velocity be- 
tween two points separated by an arbitrary distance is given by p^ : 

This leads to an energy spectrum 

2 

E{k)^'^{l~MkR)), (165) 
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which reduces to Novikov's result E{k) = ny'^/'iTTk for k 



100 . We 



observe that the velocity correlation function ( |164[ ) diverges logarithmically as 
ri —* 0. This is consistent with the divergence of the variance of the velocity 
distribution (117)-(119). Therefore, it is more proper to characterize the spatial 
correlations of the velocity by the function = (VoVi/|Vo|) which remains 

finite as ri — > 0. Its evaluation is more complex and leads to the result 



where 

is the incomplete /^-function and 



1/2 ^+cx; 



^l(^) 



— I dz. 



-HPdt, 



|Vo| 



717 

"Te" 



In TV 



1/2 



7rn IniV 



1/2 



Similarly, the correlation function K{ri) = (VoVi/Vff) is given by 



K{s) = 1 



IniV 



^ ^ ^ 4s2 ( 4s2 



where Ei{z) denotes the exponential integral 



+ 00 



Ei{z) = / dt 



(166) 



(167) 



(168) 



(169) 



(170) 



Other results concerning the characterization of the spatial velocity correlations 
can be found in lEsl. 



3.6 Statistics of fluctuations of the gravitational field 

The stochastic approach developed in the preceding sections was inspired by the 
famous work of Chandrasekhar & von Neumann |^,|^,||j24H2| concerning the 



fluctuations of the gravitational field created by a random distribution of stars. 
Let us consider a collection of N stars with mass m randomly distributed in a 
sphere of radius R with a uniform density n in average. The force by unit of 
mass created at the center O of the domain is 

N „ 

F^Y.^., = (171) 

1=1 » 

As before, the problem consists in determining the distribution of a sum of 
random variables. In the present case, the variance of the force created by one 
star diverges algebraically 

C^')-:^ K^y^-r^dr^l (172) 



r 
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so that the distribution of the total force is a Levy law. This particular Levy law 
is known as the HoUsmark distribution since it was first determined by Holtsmark 
in the context of the electric field created by a gas of simple ions |Q. In the 
thermodynamic limit N,R +00 with n — finite, the distribution of F 

can be expressed by the Fourier transform 

1 /'+00 

W(F) = ^^ ksm(kF)e-''''''^dk, (173) 



^ (2^G)3/2nTOi/2. (174) 



with 

a — 

15 

It has the asymptotic behavior 

W{F) - ^63/2^1/2^^-9/2 _^ ^^^^ (^75^ 

which can be shown to coincide with the distribution due to the nearest neighbor. 
The typical force due to the nearest neighbor Fn.n ~ Gm/d"^ ~ Gmn^^^ is 
precisely of the same order as the average value of the force due to all the stars 

(^^>=4rQ)(^)'''GW/3, (176) 

determined from Eq. ( |173D . This shows that only stars close to the star under 
consideration determine the fluctuations of the gravitational field. In fact, it is 
possible to show that the "effective" force created by a star at distance r from 
the star under consideration is ||l|: 

Gm 1 

^e// = — 1 + ,2/^2' (177) 

where yl ~ d is of the order of the interparticle distance. For weak separations, 
one has F^fj — > Gm/r^ but for large separations r ^ d, the effect of individual 
stars compensate each other and the resulting force is reduced by a factor (r/d)^. 
This corresponds to a shielding of the interaction, in a statistical sense, due to 
collective effects. 

Chandrasekhar & von Neumann have used this stochastic model to determine 
the speed of fluctuations T{F) the diffusion coefficient of stars (a calculation 
completed by Kandrup ) and the spatial |^,^ and temporal |^ correlations 
of the gravitational field etc... There is a complete parallel with the results 
obtained by Chavanis & Sire |47y48| , ^ for point vortices and this is another 
manifestation of the deep formal analogy between the two systems. 



4 Relaxation of a point vortex in a thermal bath 
4.1 Analogy with Brov^rnian motion 

We would like now to develop a kinetic theory of point vortices in order to 
describe their relaxation towards equilibrium. It has to be noted that an equi- 
librium state can be achieved in very different ways so that the kinetic theory of 
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vortices is not unique and depends on the situation contemplated (see Sec. 
In order to develop an intuition on the problem, we first propose a naive kinetic 
theory based on an analogy with Brownian theory |^,^ . The starting point of 
this analogy is to realize that the velocity of a point vortex can be decomposed 
in two terms: a smoothly varying function of position and time (V)(r,t) and 
a function V{t) taking into account the "granularity" of the system. The total 
velocity of a point vortex can therefore be written: 

V = (V)(r,i)+V(t). (178) 

The velocity (V) (r, t) reflects the influence of the system as a whole and is 
generated by the mean vorticity (uj) (r, t) according to the Biot & Savart formula 

(V)(r,t) = -i-z X I ^^-I_(^u;)ir',t)d'r'. (179) 

The fluctuation V{t) arises from the difference between the exact distribution 
inexact (r,i) of point vortices given by Eq. ( |T3| ) and their "smoothed-out" distri- 
bution {uj){r,t) = 717. Indeed, if we consider a surface element a, the number 
of point vortices actually present in this area will fluctuate around the mean 
number na. These fluctuations will be governed by a Poisson distribution with 
variance na. On account of these fluctuations, the velocity of a vortex will depart 
from its mean field value (V). The velocity fluctuation V, of order 7/d (where 
the inter- vortex distance), is much smaller than the average velocity 
{V), of order n"fR (where R is the domain size), but this term has a cumulative 
effect which gives rise to a process of diffusion. It makes sense therefore to in- 
troduce a stochastic description of the vortex motion such as that for colloidal 



suspensions in a liquid |81 2^ or stars in globular clusters |2^. However, con- 
trary to the ideal Brownian motion, point vortex systems have relatively long 
correlation times so that V{t) is not a white noise. This makes the study much 
more complicated than usual and the technical developments of Sec. ^ are re- 
quired. However, in order to gain some physical insights into the problem, we 
shall ignore this difficulty for the moment and describe the system by traditional 

stochastic processes. 

According to Eq. (|l78| ), we would naively expect that the evolution of the 
density probability P(r, t) would be governed by a diffusion equation of the form 

^ + (V>VP = DAP, (180) 

where D is the diffusion coefficient. This would in fact be the case for a passive 
particle having no retroaction on the vortices or when the distribution of vortices 
is uniform in average like in Sec. 0. However, this diffusion equation cannot be 



valid when the system is inhomogeneous. Indeed, Eq. (18C) does not converge 
towards the Boltzmann distribution (^ ) when t — s- -l-oo. It seems therefore that 
a term is missing to act against the diffusion. 

This problem is similar to the one encountered by Chandrasekhar in his 
stochastic approach of stellar dynamics Chandrasekhar solved the prob- 
lem by introducing a dynamical friction in order to compensate for the effect 
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of diffusion. The occurence of this frictional force is a manifestation of the 
"fluctuation-dissipation" theorem in statistical mechanics. In the present con- 
text, the dynamical friction is replaced by a systematic drift of the vortices pO| . 
This drift appears when the distribution of vorticity is inhomogeneous. It can be 
understood in terms of a polarization process like in a phenomenon of induction. 



In Sec. 4.3, we shall derive the drift term directly from the Liouville equation 
by using a linear response theory. For the moment, we introduce this term by 
hands and rewrite the decomposition ( |l78| ) in the form 

V= (V) -eVV^+V(t), (181) 



where ^ is the drift coefficient. Equation (181) must be viewed as a stochastic 
equation analogous to the Langevin equation in the ordinary Brownian theory. 
The corresponding Fokker-Planck equation can be written [pof : 

dP 

— + {V)VP^V(DVP + ^PVip). (182) 
ot 

The physical interpretation of each term is straightforward. The left hand side 
(which can be written dP/dt) is an advection term due to the smooth mean 
field velocity (V). The right hand side can be written as the divergence of a 
current — VJ and is the sum of two terms: the first term is a diffusion due to the 
erratic motion of the vortices caused by the fluctuations of the velocity (see Sec. 
H) and the second term accounts for the systematic drift of the vortices due to 
the inhomogeneity of the vortex cloud. At equilibrium, the drift precisely bal- 
ances random scatterings and the Boltzmann distribution (|4^ ) is set tled . More 



precisely, the condition that the Boltzmann distribution satisfies Eq. (182) iden- 
tically requires that D and ^ be related according to the relation 

^ = DPj, (183) 

which is a generalization of the Einstein formula to the case of point vortices. It 
is remarkable that we can obtain such a general relation without, at any point, 
being required to analyze the mechanism of "collisions" . A more rigorous justi- 
fication of this relation will be given in the next sections in which the diffusion 
coefficient and the drift term are calculated explicitly. Note that the diffusion 
current can be written J — xVa where a = ln(ti;) -f f3"fip is a "generalized 
potential" which is uniform at equilibrium, see Eq. (H). Therefore, the Fokker- 



Planck equation (182) is consistent with the linear thermodynamics of Onsager 



which relates the diffusion currents to the gradients of generalized potentials. The 



Fokker-Planck equation (182) can also be obtained from a variational formula- 
tion (the so-called Maximum Entropy Production Principle): it can be shown to 
maximize the rate of entropy production 5* under the constraints brought by the 
dynamics Q. 

4.2 Diffusion coefficient of point vortices in an inhomogeneous 
medium 

In this section, we determine the value of the diffusion coefficient D which enters 



in the Fokker-Planck equation (182). The essential difference with the calculation 
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of Sec. ^ is that the vorticity distribution is now inhomogeneous so that there 
exists a shear in the system. This shear considerably modifies the expression of 
the diffusion coefficient [ p0||3^ . An explicit expression for D can be obtained 
in a "thermal bath approximation" which is valid if we are sufficiently close to 
equilibrium. To be specific, we consider a system of N point vortices at statistical 
equilibrium with an inverse temperature (3eq- These "field vortices" form the 
thermal bath. We now introduce a "test vortex" in the system an study its 
diffusion in the sea of "field vortices" . For simplicity, we restrict ourselves to 
simple equilibrium flows which are either unidirectional or axisymmetric. 

In the case of unidirectional flows in the x-direction, we define the diffusion 
coefficient in the y direction by 

D = limt^+^^^^p- with Ay^ I Vy{t')dt', (184) 

Jo 

where V(t) is given by Eq. (M). It is possible to put the diffusion coefficient in 
the form of a Kubo formula [ |36{ : 

p + CO 

D= {Vy{t)Vy{t^T))dT, (185) 

Jo 

where the quantity in bracket is the velocity auto-correlation function at different 
times. Using Eq. (|lj), we have explicitly 

p + CO P 

D = N j dT J d^riVyil^O,t)Vyil^O,t-T)P,,iyi). (186) 
For axisymmetric equilibrium flows, one has similarly 

\2\ pt 



D = limf^+oo-^^^^ with Ar= I Vr{t')dt', (187) 
2^ Jo 



and 



r+OD p 

D = nJ^ dry d2riK(t)(1^0,i)K(*-r)(1^0,t-r)Pe,(ri). (188) 

The velocity auto-correlation function is calculated in Appendix A under the 
assumption that between t and t ~ t the point vortices follow the equilibrium 
streamlines, which is valid for sufficiently strong shears. It is found that the 
auto-correlation function decreases like at large times so that the diffusion 
coefficient is well-defined. It must be stressed, however, that the decorrelation is 
slow so that the fiuctuations cannot be described by a white noise process, as 
indicated previously Q Using the results of Appendix A, the final expression for 
the diffusion coefficient in the presence of a shear can be put in the form [3C|j3^ : 

^=|^^lniV(^),„ (189) 



^ The problem is even more severe in stellar dynamics where the temporal correlation 
function of the gravitational force decreases like t~^, responsible for logarithmic 
divergences in the diffusion coefficient |88]. 
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where |£'(r)| is the local shear of the flow: \S\ = —d/dy{V)eq{y) for unidhec- 
tional flows and \IJ\ — rd/ dr{{V) eq(r) / r) for axisymmetric flows. In Sec. 5 



we 

shall calculate the diffusion coefficient in another manner and extend formula 



(189) to more general equilibrium flows with |£'(r)| = 2-\/— det(Z') where E is 
the stress tensor. The expression (189) for the diffusion coefhcient can be written 
in the general form (see Sec. 0): 



(190) 



with a correlation time r = 27r/|Z'(r)|. Physically, it corresponds to the time 
needed by two vortices with relative velocity Sd (where d is the inter-vortex 
distance) to be stretched by the shear on a distance ~ d. Of course, when there 
is no shear, our approximations break down and the diffusion coefficient is given 
by Eq. (154) which involves a correlation time r ~ 1/ {uj){r)\^lnN determined 
by the dispersion of the vortices |^ . These two results correspond to a limit of 
"strong shear" and "weak shear" respectively. Clearly, a general formula for D 
should take into account simultaneously the effect of the shear and the dispersion 
of the vortices. 



4.3 Systematic drift experienced by a point vortex in an 
inhomogeneous medium: linear response theory 



In Sec. 4.1, we have argued on the basis of general considerations that when 
the vorticity distribution is inhomogeneous, a point vortex should experience a 
systematic drift in addition to its diffusive motion. In this section, we justify this 
drift by a linear response theory, starting directly from the Liouville equation 
Consider a collection of N point vortices at statistical equilibrium. In the 
mean-field limit, the N-particle distribution function /Lteg({rfe}) can be approxi- 
mated by a product of N one-particle distribution functions Peq, each of which 
at equilibrium with the same inverse temperature f3eq- 

N N 

f^eqiUk}) = n ^-^(rfe) = n Ae'P'^''^'^(^^\ (191) 
fc=i fc=i 

This distribution is stationary, in a statistical sense, since it corresponds to a 
maximum entropy state. The introduction of an additional point vortex, the 
"test vortex" , will modify this equilibrium state. The distribution function of 
the field vortices becomes time dependant and can be written in the form 

fi{{rk}, t) = fieqiUk}) + M'({rfc}, t), (192) 

where the perturbation /x'({rfc},i) reflects the influence of the test vortex on 
its neighbors, just like in a polarization process. The A''-particle distribution 
function /i({rfe},i) satisfies the Liouville equation 



N 



dt 
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Substituting Eq. ( |192| ) in Eq. ( |193D , we obtain the evolution equation of the 
perturbation /i': 



N 



dipeq , 



(194) 



where C = J2iLi ^ Liouville operator and V.^ ~ 'l2j^i ~^ i)+y{0 

i) denotes the total velocity of vortex i. This equation can be solved formally 
with the Greenian 



G{t,t') = expj-^ /:(r)dr|. 



(195) 



If i = is the time at which the test vortex is introduced in the system, we have 
fj,'{t^O) = 0. One then finds that 



M'(t)=/?e,7 f drGit,t-r)J2y.^{ 
Jo 



r,)fleqi{rk}). 



(196) 



The average velocity of the test vortex is expressed in term of the distribution 
function /i of the field vortices by 



N 

V) = / n^'r,V/i({rfe},i), 

k=i 



(197) 



where V = J2f=i V(z ^ 0). Inserting Eqs. (|192|) and (|196| ) into Eq. ( |197D , one 
obtains 

(V)= / l[d\kVfieq{{rk}) + f3eql / n^'^'^V 

k=l •' k=l 

X rdrG(i,t-r)f]y,''^(r,)A*e,({r,}), (198) 

with summation over repeated greek indices. The two terms arising in this ex- 
pression have a clear physical meaning. The first term is the mean field veloc- 
ity (V)eq — — z X VV'eg(r) Created by the unperturbed distribution function 
Meij({rfe}). The second term, arising from the perturbation /i', corresponds to 
the response of the system to the polarization induced by the test vortex. Be- 
cause of this back reaction, the test vortex will experience a systematic drift 
(y) drift = (V) — (y)eq- Explicating the action of the Greenian ( |195| ), we obtain 



(V)drift = Peql 
+ V''i0^i,t-T) 



[ ^d'r,5^V(^^0,^) / dTj2 T.^'^ij 

k=l i=l i=l '-j^i 



d^peq 



^(r,(i-r))Me,({rfe(t-r)}), 



(199) 
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where ri(t — r) is the position at time t — r of the point vortex i located at 
Ti (t) = at time t. This is obtained by solving the Kirchhoff-Hamilton equations 
of motion dvi/dt = between t and t — t. 



The exact expression of the drift (19£) is completely inextricable in the gen- 



eral case. In order to enlighten its physical content, we have to make some 
approximations. We shall consider in the evaluation of the time integral that the 
point vortices are advected by the equilibrium mean- field velocity (y)eq- This 
is reasonable in a first approximation because, when N — s- cx), the typical ve- 
locity fluctuations V, of order ^ /d ^ ^N^^^ /R are much smaller than the mean 
field velocity (V)eg, of order N-f/R. Of course, this approximation breaks up at 
scales smaller than 5 ~ R/N when the velocity fluctuations become comparable 
to the average velocity. In that case, we cannot ignore the details of the discrete 
vortex interactions anymore and a specific treatment is necessary. For simplic- 
ity, we shall introduce a small-scale cut-off and replace the exact Greenian G 
by a smoother Greenian {G)eq constructed with the averaged Liouville operator 
{C)eq = J2iLi{^^)eq'^- t^is approximation the correlations involving two 
different vortex pairs vanish and we obtain 

{V^)dr^ft^^3eql / l[dhk / dT^y^(z^O,t) 

fc=i -^o t=l 



xr^(0 -> - T)^{Mt - T))Y[P,,{rk), (200) 



N 



dr: 



i—k 



where we have used P'^'^{rk{t — t)) = P'^'^{rk{t)) since Peg — f{^eq) is constant 
along a streamline. Since the vortices are identical, we have equivalently 



{V^)dr^^t = Npeql / d'v^ / dTV^'{l 0, i) 







xl/-(0->l,t-r)^(ri(i-r))P,,(ri). (201) 
In the case of a unidirectional equilibrium flow in the x direction, the drift 



velocity (201) can be written 



{Vy)dr^ft^N^3eql J d^l J dTVy{1^0,t)Vy{Q^l,t-T)^{yi)Peq{yi), 

(202) 

where we have used yi {t — t) — yi (t) = yi and the time integration has been 
extended to +oo. Since the space integration diverges as ri — > Tq, we can make 
a local approximation dy^{yi) ~ dyip{y), neglect the contribution of images in 
the velocity Kernel and use Vy{0 1) = —Vy{l 0). The local approximation 
reflects the strong influence of the nearest neighbor and is only marginally valid 
as discussed in Sec. |3[ The expression of the drift can then be written 

{Vy)dr^ft = -PeqlD^iy), (203) 
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where D is the diffusion coefficient given by formula (186). The same type of 
relation is obtained for an axisymmetric equilibrium flow. More generaly, one 
can write the drift as | |30| , p6[ : 

{V) drift = -fieqlDW^^q, (204) 



where D is given by Eq. ( [1891) . We have thus derived, from a kinetic theory, an 
Einstein relation ^ = D^Peq for the point vortex gas. This relation is of great 
conceptual importance and could be checked by direct numerical simulations of 
point vortex dynamics in the thermal bath approach. 

The direction of the drift has also important physical implications. First, we 
note that the drift is always normal to the mean field velocity (V)eq — —zx'V^'eq- 
Since Peg = fitpeq) for an equilibrium flow, this implies that the drift is directed 
along the vorticity gradient. This is clear at first sight if we write the expression 



of the drift in the form (V)drift — DV ln{u>)eq, using Eqs. (|204| ) and (|44|). Since 
D > 0, this formula indicates that the test vortex always ascends the vorticity 
gradient. In fact, we have assumed in the previous discussion that all the vortices 
have the same circulation. It is straightforward to generalize our calculations for 
a test vortex with negative circulation evolving in a bath of vortices with positive 
circulations. Similar results are obtained except that now the test vortex descends 
the vorticity gradient. 

If we now take into account simultaneously the drift and the diffusion of the 
test vortex, we can argue that the evolution of the density probability P(r, t) is 
governed by the Fokker-Planck equation 

dP 

— + (V)e,VF = V(Z?(VP + PeqlP'^^eq)), (205) 

where we recall that ij^eq is the stationary streamfunction generated by the equi- 
librium vorticity {'jj)eq via the Poisson equation (^. In the case of a unidirec- 
tional flow, this Fokker-Planck equation can be transformed into a Schrodinger 



equation (with imaginary time) which can be solved analytically |116| 



4.4 A relaxation equation for point vortices 

If we are sufficiently close to equilibrium, we can try to apply our previous results 
to all vortices in the system, eliminating the somewhat arbitrary distinction 
between test and field vortices. We propose to describe the relaxation of a cloud 
of point vortices towards statistical equilibrium by the following set of equations 

dP 

— + {Y)VP ^V{D[VP + fi-iPVi')), (206) 
at 

= ^N-fP, (207) 



consisting of the Fokker-Planck equation (205) coupled to the Poisson equation 
(0). This model is expected to be valid only close to equilibrium so that the 
inverse temperature P has a clear physical interpretation. In this model, a point 
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vortex is assunied to undergo a diffusion process due to the fluctuations of the 
velocity and a systematic drift (V) drift — —(i^DVijj due to the inhomogeneity of 
the vortex cloud. At negative temperatures, the drift is directed inward and the 
vortices tend to duster while at positive temperatures, the drift is directed out- 
ward and the vortices tend to repell each other and accumulate on the boundary. 
For /? = 0, there is no drift and Eq. (p06|) reduces to the pure diffusion equation 



(180); in that case, the vorticity distribution is uniform in aver age. These results 



are consistent with the thermodynaniical approach of Onsager 1 101 1 but the drift 
provides a physical mechanism to understand the clustering of point vortices at 
negative temperatures. It can also be noted that Eq. ( ^06[) is formally similar to 
the Smoluchowski equation describing the relaxation of colloidal suspensions in 



an external gravitational field |108]. In the present context, however, the field ip 
is not fixed but is generated by the distribution of particles itself via the Pois- 
son equation. The resulting Smoluchowski-Poisson system has been studied in 



6yilq| for various space dimensions. 



The system of equations ( |20q ) and (207) is also similar to the model intro- 



duced by Chandrasekhar in his stochastic description of stellar dynamics: 

= 47rG J fd^v. (209) 

In that model, a star undergoes a diffusion process (in velocity space) due to 
the fluctuations of the gravitational force and a dynamical friction (F) friction = 
—DPmv resulting from close encounters. Fundamentally, this friction is due 
to the inhomogeneity of the velocity distribution. The coefficient of dynamical 
friction is given by an Einstein formula ^ = D[3m in which the velocity dispersion 
1 //? of the stars enters explicitly. 



The morphological similarity of the two models (p0q)-(|20^) and (|08|)-(|0g) 
is striking although the physical content of these equations is, of course, very 
different. In this analogy, we see that the systematic drift of the vortices is the 
counterpart of the dynamical friction of stars. In fact, Chandrasekhar's dynam- 
ical friction can be derived from a linear response theory [fz^ exactly like we 
have derived the systematic drift of a vortex. In addition, both terms can be un- 
derstood physically as a result of a polarization process (see discussion in |3^). 
This is another mark of the formal analogy between point vortices and stars. 



5 Kinetic theory of point vortices 

The previous relaxation equations are only valid for a test vortex evolving in a 
bath of field vortices, or for a cloud of point vortices close to statistical equi- 
librium. We would like now to relax this "thermal bath approximation" and 
describe more general situations which do not explicitly rely on the existence of 
a well-defined temperature or equilibrium state. Hence, we would like to develop 
a more complete kinetic theory of point vortices [p6[ . 
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5.1 The Liouville equation 

Let us consider a collection of + 1 point vortices with identical circulation 7. 
Let /i(r,ri, ...,rjv,t) denote the + 1 particle distribution of the system, i.e. 
/u(r, ri, r^r, t)(PrcPri...(PrN represents the probability that point vortex be 
in the cell (r, r + dr), point vortex 1 in the cell (ri, ri + dri)... and point vortex 
N in the cell (rA?, ta? -|-drAr) at time t. The (A^ + l)-particle distribution function 
IJ,{t) satisfies the Liouville equation 

where is the velocity of vortex i produced by the other vortices according 
to Eq. (p^. We also introduce the one- and iV-particle distribution functions 
defined by 

N 



P(r,t)= / fi{{rk},t)f[d^rk, (211) 

k=l 

Ai.j,.(ri,...,rjv,i) = J fi{{rk},t)d^r. (212) 
We write the distribution function fi in the suggestive form 

^(r,ri, ...,rN,t) = P{r,t)fisys{ri, ■■■,rN,t) + ^/(r, ri, ■■■,rN,t), (213) 
where the quantity fij refiects the effect of correlations between vortices. 



The Liouville equation (210) provides the correct starting point for the analy- 
sis of the dynamics of our vortex system. However, when iV is large, this equation 
contains much more information than one can interpret. Consequently, what one 
would like to do is to describe the system in some average sense by a one-particle 
distribution function. 



5.2 The projection operator formalism 

Our first objective is to derive some exact kinetic equations satisfied by P{r, t) 
and /isj^s(ri, rjv, t). This can be achieved by using the projection operator 



formalism developed by Willis & Picard |127]. This formalism was also used by 
Kandrup | }73[ in the context of stellar dynamics to derive a generalized Landau 
equation describing the time evolution of the distribution function of stars in 
an inhomogeneous medium. We shall just recall the main steps of the theory. 



More details can be found in the original paper of Willis & Picard [127| and in 
Kandrup ||73|]. T o have similar notations, we set x = {r} and y = {ri, ...,rjv}. 
Then, Eq. ( ^13| ) can be put in the form 

Kx, y, t) = hr{x, y, t) + mix, y, t), (214) 



with 



m{x,y,t) = f{x,t)g{y,t), 



(215) 
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where we have written f{x,t) = P{r,t) and g{y,t) = /isj,s(ri, r^r, t). The 
Liouville equation is also cast in the form 



'di 



— = -iLfi^ ~i{Lo + Lsys + L')fi, (216) 



where Lq and L^ys act respectively only on the variables x and y, whereas the 
interaction Liouvillian L' acts upon both x and y (the complex number i is here 
purely formal and has been introduced only to have the same notations as in 



Following Willis & Picard |127|, we introduce the time-dependant projection 
operator 



P{x, y, t) = g{y, t) J dy + f{x, t) J dx - f{x, t)g{y, t) J dx J dy. (217) 

We can easily check that 

P(x, y, t)n{x, y, t) = fiB^x, y, t), (218) 

[1 - P{x, y, t)]^l{x, y, t) = ^/(x, y, t). (219) 

We also verify that P is a projection in the sense that P^(i) — P{t). Applying 
P and 1 — P on the Liouville equation (pl6|), we obtain the coupled equations 



^t^J■R{x, y, t) = -iPLfiR - iPLfii, (220) 

and 

dtfiiix, y, t) = -^(l - P)LfiH - i{l - P)Lfij. (221) 

These equations describe the separation between a "macrodynamics" and a "sub- 
dynamics" . 

Introducing the Greenian 

g{t, t') = exp|-i ^ dt"[l - P{t")]L^ , (222) 



we can immediately write down a formal solution of Eq. ( p2l| ) , namely 

^lIix, y,t) = - f dt'Git, t')i[l - P{t')]LfiRix, y, t'), (223) 
Jo 

where we have assumed that, initially, the particles are uncorrelated so that 



yU/(x, y, 0) = 0. Substituting for fJ.i{x, y, t) from Eq. (223) in Eq. (220), we obtain 

dtm{x,y,t) = -iPL^iR- f dt'P{t)Lg{t,t')[l-P{t')]LfiR{x,y,t'). (224) 

Jo 
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The integration over y will yield an equation describing the evolution of /. Using 
some mathematical properties of the projection operator (217), the final result 
can be put in the nice symmetrical form given by Willis & Picard |127|, 

dtf{x, t) + iLof + i{L')sysf = -J dt' J dyAtL'Git, t')At,L'g{y, t')f{x, t'), 

(225) 

(226) 

(227) 
(228) 



where the notations stand for 

{L')sys = J dy'L'{x,y')g{y',t), 

{L')o = I dx'L\x\y)f{x\t), 

AtL' = L' -{L'),ys-{L')o. 
Similarly, after integrating over x we find the equation satisfied by g, 



dtg{y,t) + iL,y,g ~ i{L')ig ^ - / dt' / dxAtL'g{t,t')ArL'g{y,t')f{x,t'). 

(229) 



5.3 Application to the point vortex system 

The previous theory is completely general and we now consider its application to 
a system of point vortices [ p6[ . Let us first rewrite the Liouville equation ( 210 ) in 
a form that separates the contribution of the test vortex from the contribution 
of the field vortices: 

| + i:v(.-o,|;+x;v(o^,-)|ii 



4 = 1 

N N 



dr 



E E v(,^^)l^^o. 



j=i 



(230) 



Applying the general theory of Willis & Picard |127], we obtain the following 
kinetic equation for the one-particle distribution function of a vortex system : 



dP 
'dt 



(V 



dP 
dr 



a 



t „ N N N 

/ ^-/n^^r^-EE^^( 



0) 



xg(i, t - r) ( V''(j -> 0)— + V"^(0 ^ j)— )P(r, t - T)fisy4{rk}, t - r). 



dr''- 



(231) 

where the Greek indices refer to the components of V in a fixed system of co- 
ordinates and V{i — > 0) = V(i 0) — (V(i 0)) denotes the velocity fiuctu- 
ation. We can note that equation ( |23l| ) already shares some analogies with the 
Fokker-Planck equation of Sec. ^ Indeed, the first term on the right hand side 
corresponds to a diffusion and the second term to a drift. For a passive particle 
V'^(0 j) = and the drift cancels out, as expected. 
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5.4 The factorization hypothesis 

If the vortices are initially decorrelated then, for sufficiently short times, they will 
remain decorrelated. This means that the {N + l)-particle distribution function 
can be factorized in a product of {N + 1) one-particle distribution functions 

N 

^i{r,ru-,rN,t)^l[P{rk,t). (232) 

k=0 



If we integrate the Liouville equation (21C) on the positions of the N vortices 
1, and use the factorization ( ^32|) , we directly obtain |36| : 

dP 

— + (V>VP = 0. (233) 

Therefore, for sufficiently short times, the average vorticity (w) satisfies the 2D 
Euler equation (called the Vlasov equation in other circumstances). However, 
at later times, the distribution function /i differs from the pure product ( ^32| ) 
and the Euler equation does not provide a go od ap proximation anymore. In Sec. 
^.3| we have determined an exact equation (231) satisfied by the one-particle 
distribution function which is valid at any time. This equation is not closed, 
however, since it involves the A'^- vortex distribution function ^sys- We shall close 
the system by assuming that pLgys can be approximated by a product of N one- 
particle distribution functions in the form 

JV 

t)^\{P{vk,t). (234) 
fc=i 



0) 



Inserting Eq. ( ^34[ ) in Eq. ( |231| ) , we obtain 

•' k = l ■^'^ 4 = 1 7 = 1 



J k—1 



N 



xg{t, t-T)[ V^U ^Q)±-+V%Q^3)^] P{v, t - r) n P{^k,t - t). 



(235) 



If we assume that between t and t ~ t the trajectories of the particles are de- 
termined from the smooth velocity field created by the vorticity distribution 
(oj) — N'^P{r,t), the foregoing equation simplifies in 

X ^ 0)P,^ + l'-(0 ^ 1)-P|^}__^. (236) 
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where P — P{r,t) and Pi = P{ri,t). Eq. (236) is a non Markovian integrod- 
ifferential equation since the probabihty density P(r, t) in r at time t depends 
on the value of the whole distribution of probability P(ri, t — r) at earlier times 
through an integration on ri and r. Equation ( |236| ) is therefore non local in space 
and time. It can be shown |36 that this kinetic equation rigorously conserves 
angular momentum in a circular domain and linear impulse in a channel (or in 
an infinite domain). However, under this form, it has not been possible to prove 
the conservation of energy and the H-theorem for the Boltzmann entropy (^l|). 



5.5 The case of short decorrelation times 

If we assume that the decorrelation time r is short (which does not need to be 
the case) and implement a strong Markov approximation, we obtain 

x|f^(1 ^ 0)Fi^ + V^{0 ^ 1)P|^}. (237) 

In the case of an infinite domain V(0 ^ 1) = — V(0 1) and we have the 
further simplification 



where 



K'^-iO = ^ = g4 (239) 



and ^ = ri — r. To arrive at Eq. ( |239| ) we have explicitly used the form of the 
Kernel (^4|), and to get the second equality we have used the f act that we are in 
two dimensions. Note that the symmetrical form of Eq. ( 238 ) is reminiscent of 
the Landau equation introduced in plasma physics and in stellar dynamics (see, 
e.g., II): 



K^^'^ 2TTNG^m^ ln( ) "''^"'^ ^ """'^ , (241) 



with 



and u = V — vi, / = /(v, t), fi — /(vi, <). In this analogy, the position r of the 
vortices plays the role of the velocity v of the electric charges or stars and the 
spatial distribution P(r, t) the role of the velocity distribution /(v, t). Therefore, 
we can directly infer the conservation of linear impulse = J{uj)r(Pr and 
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angular momentum L = J{uj)r^(Pr which play respectively the role of impulse 
P = J fv(Pv and kinetic energy K = J f^d'^-v in plasma physics. We can also 
prove a 77-theorem for the Boltzmann entropy ( p| ) exactly like in the case of 



the Landau equation. Finally, the solutions of Eq. (238) converge towards the 
Gaussian vortex (the equivalent of the Maxwellian distribution in plasma physics 
with r in place of v): 

P(r) = Ae-5"7(r-ro)^ (242) 

which is the maximum entropy state at fixed circulation, angular momentum 
and impulse. It is in general different from the Boltzmann distribution ( ^^ 
with the relative streamfunction ip' — ij: + — Uj^r, except in the particular 
limit /3 ^ 0, 12 — > +oo with fixed a — (3tl/2 corresponding to the statistical 
equilibrium (|8^) . This clearly indicates that Eq. ( |23§| ) does not conserve energy. 

5.6 A generalized kinetic equation 

Now, if we account properly for memory effects in Eq. ( |236| ) , we can obtain 
a generalized kinetic equation which guaranties the conservation of energy (in 
addition to the other constraints) and is therefore more satisfactory. If the distri- 
bution of vortices is axisymmetric, it is possible to calculate the memory function 
appearing in Eq. ( |236| ) explicitly if we assume that the correlation time is smaller 
than the typical time on which the average vorticity changes appreciably | |3^ . 
In this approximation, the point vortices follow, between t and t — r, circular 
trajectories with angular velocity f2{r,t) = {Vg){r,t)/r and Eq. (|236| ) simplifies 
in (see Appendix B): 

at Ar or Jq L V^>/ J ''i '^^i 

(243) 



where fi — n{r, t), Qi — Q{ri,t) and r> (resp. r<) is the biggest (resp. smallest) 
of r and ri . The angular velocity is related to the vorticity by 

{u) = ~inr^). (244) 
r or 

We can propose an approximation of the general kinetic equation (|236| ) which 
encompasses the axisymmetric form previously derived. Memory effects are not 
neglected, unlike in Eq. ( ^38[ ), but they are simplified in a way which preserves 
all the conservation laws of the system (as discussed below). We propose the 
generalized kinetic equation pq] : 



■ (V)VP - ' 



dt ' ' 8 ^r^' 

with 



/d^r,i^-J(ev)(p,^-P^), (245) 



^'^^(0 = -^ = (246) 
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and ^=ri — r,v = (V) (r i , <) — (V) (r , t) . In the thermodynamic hmit N +00 
and 7 ^ 1/iV (see Sec. 2.7), the kinetic equation (245) reduces to the Vlasov 
equation ( ^331 ). However, in practice, N is always finite and the correlations 
between point vortices must be taken into account. The "collision term" in Eq. 
( ^45D gives the first order correction 0(1/ N) to the Vlasov limit. 

From Eqs. ( ^43|) and (245), it is clear that the relaxation of the point vortices 
is due to a phenomenon of resonance. Only the points ri satisfying the condition 
^•v = with ri 7^ r contribute to the diffusion current in r. In the axisymmetric 
case, this condition of resonance reduces to J7(ri) = f2{r), which supposes that 
the angular velocity profile is non monotonic. The occurence of a 5-function in 
Eqs. ( ^43| ) and ( |245| ) is the main difference with the Landau equation (240). In 
the present context, it ensures the conservation of the "potential" energy of the 
vortices E = ^ J ui'ip(Pr which has no counterpart in Landau's theory applying 
to spatially uniform plasmas. It can also be noted that, contrary to the Landau 
equation, the kinetic equations (23S) and (245) do not suffer the well-known 
logarithmic divergence appearing in the context of Coulombian plasmas and 
stellar systems (see, e.g., ||7^). This is due essentially to the lower dimension of 
space (D = 2 instead of D = 3) and to the different nature of the interactions. 



5.7 Conservation laws and H-theorem 

We now derive the conservation laws and the H-theorem satisfied by Eq. ( |245| ). 
The conservation of the circulation is straightforward since the right hand side 
of Eq. ( |245| ) can be written as the divergence of a current. To prove the conser- 
vation of angular momentum, we take the time derivative of L = Nj J Pr^cPr, 
substitute for Eq. ( |245| ) , permut the dummy variables r and ri and add the 
resulting expressions. This yields 

i = ^/ .»rf,A-."^..(ev,(p,|^ - P|^i). ,247) 

From Eq. ( p46| ) , we immediately verify that 

K^"'^'' = 0, (248) 

which proves the conservation of angular momentum. We can prove the conser- 
vation of linear impulse in a similar manner | |3^ . For the conservation of energy, 
we start from Eq. ( p9| ) and follow the same procedure. This yields 



^ = '^'rd^riif-<5(ev)(^A— -P^). (249) 



Considering the form of the tensor (246), we have 



K^'vl = p(ev). (250) 
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When substituted in Eq. (249), we see that the occurence of the (5-function in 
the kinetic equation implies E = 0. Finally, for the rate of entropy production 
we have, according to Eqs. 



l|) and ( p45D 
1 



(251) 



Permutting the dummy variables r and ri and adding the resulting expression 
to Eq. (|251|), we obtain 



5 = 



16 



1 



PPi 



dP 



dPi 



9< 



p^dp_pm 



(252) 

Now, for any vector, A^^K^^''A'' = (A •Ci)^/C^ > 0. This proves a H-theorem 



{S > 0) for the kinetic equation (245). It should be emphasized that the con- 
servation laws and the H-theorem result essentially from the symmetry of the 
kinetic equation. This is satisfying from a physical point of view. 
It is also easy to show that the Boltzmann distribution 



P 



is a stationary solution of Eq. ( |245[ ) . Noting that 



(253) 



(254) 



we have successively 

f BP BP \ 

\^'B?^ " ^~B^J ^ /57mif'^^« + - /37m|^(ev), (255) 

wher e we have used Eqs. ( |24§| ) and ( |25C| ). When Eq. (|255| ) is substituted in 
Eq. (245), we find that the right hand side cancels out due to the (5-function. 
The advecti ve te rm is also zero since P = f{tp '). T herefore, the Boltzmann 
distribution (253) is a stationary solution of Eq. ( ^45[ ) . Note, however, that this 
is not the only solution, unlike for ordinary kinetic equations. Any stationary 
solution of the Euler equation satisfying in addition ^.v 7^ for any couple of 
points r, ri (with r ^ ri) is a solution of Eq. ( ^45| ). Physically, this implies that 
the system needs sufficiently strong resonances to relax towards the maximum 
entropy state. If this condition is not realized, the system can remain frozen in 
a sort of "metastable" equilibrium state. Further evolution of the system will 
require non trivial correlations between point vortices which are not taken into 
account in the present theory. 



5.8 The thermal bath approximation 

A direct connexion between the generalized kinetic equation ( |245| ) and the 
Fokker-Planck equation of Sec. ^ can be found. Introducing a diffusion tensor 

f j2„ r^u•^,r^,.^T, (256) 



Jd\,K^^''S{^.^)Pi, 
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and a drift term 



7J^ = 



Nr 



we can rewrite Eq. (245) in the more illuminating form 



dP 



(257) 



(258) 



similar to a general Fokker- Planck equation. Note, however, that Eq. ( ^58[ ) is an 
integrodifferential equation since the density probability P(r, t) in r at time t 
depends on the value of the whole distribution of probability P(ri, t) at the same 
time by an integration over ri . By contrast, the Fokker-Planck equation ( 182 ) is a 
differential equation. The usual way to transform an integrodifferential equation 
into a differential equation is to make a guess for the function -P(ri) appearing 
under the integral sign and refine the guess by successive iterations. In practice 
we simply make one sensible guess. Therefore, if we are close to equilibrium, 
it seems natural to replace the function Pi appearing in the integrals by the 
Boltzmann distribution 

P(ri) = Ae^'^'^'^f''!). (259) 

This corresponds to the "thermal bath approximation" of Sec. |4[ the vortices 
have not yet relaxed completely, but when we focus on the relaxation of a given 
point vortex (described by P) we can consider, in a first approximation, that the 
rest of the system (described by Pi) is at equilibrium. Within this approximation, 
the diffusion coefficient and the drift simplify in 



P(r, t) J K^"'6{t^r)d% 



(260) 
(261) 



where we have made the local approximation. If we assume that the correlation 
time is short, i.e. if we replace £,^S{^.v) by t/tt^ (compare Eqs. (245) and ( |238| )), 
we obtain 

= PjDV^, (262) 



lOTT 



and Eq. (258) reduces to the Fokker-Planck equation 

dP 
'dt 



(V>VP = \7iD{\7P + (3-fPV^Jj)). 



(263) 



(264) 



This approximation is, however, not very satisfactory since the decorrelation 
time r appears as a free parameter. In fact, the decorrelation time can be de- 
termined self-cons isten tly from the above formulae by evaluating properly the 
5- function in Eq. (261). Expanding the velocity difference v = (Vi) — (V) in a 
Taylor series in ^ = ri — r, we obtain to first order in the expansion 



(265) 
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where 

E^^^U'4r^mi\ (266) 

is the stress tensor. It has the property of symmetry Z''"' — U'^^. Since th e flow 
is divergenceless, we also have JJ^'^ + JJ^y = 0. In terms of the stress tensor (266), 
the diffusion tensor ( |26l| ) can be rewritten 

D'"' = — ^ J ^ ^ ^ S{S^"'eC)d% (267) 

This integral can be performed easily by working in a basis {^[, ^2) which the 
tensor 17^" is anti-diagonal Then, 

D^" = j ^:^^5mm'iQ d^'^, (268) 

where we have set |£'(r)| — 2-^/— det(Z'). Clearly, this quantity is invariant by a 
change of reference frame and it measures the local shear of the flow. It is easy 
to check that the diffusion is isotropic and that = D5^^'^ with 



Setting = ^cos6' and ^2 — ^sin0 where ^ = ^' = |ri — r|, we obtain 

L>=— ^— — / d9sm^ed{^^cos9sme), (270) 

8 l-^lr)! Jo Jo 



or, equivalently, 

,,2 



As explained previously, we regularize the logarithmic divergence by introducing 
appropriate cut-offs at small and large scales. With the change of variables t = 
cos 9, we finally obtain 



which establishes Eq. ( |189D in the general case. 



5.9 The collisional relEixation time 

We can deduce from this kinetic theory the "collisional" rel axati on time of the 
point vortex gas. Considering the Fokker-Planck equation (264), it is easy to 
check that, for t +00, the distribution function P(r, t) will converge towards 
the Boltzmann distribution (E3). The relaxation time corresponds typically to 
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the time needed by a vortex to diffuse over a distance R, the system size. There- 
fore, treiax ~ / D , with Z) ~ 7 In TV according to Eq. ( p72\) . Using F = N"f and 
introducing the dynamical time tjj ^ i^^)^^ ^ /T, we obtain the estimate 

N 

treiax ^ 1 77^0, (273) 

In A* 

as in the case of cohisional stehar systems . Since the statistical description is 
expected to yield relevant results for large N, we conclude that the "coUisional" 
relaxation of point vortices towards the Boltzmann distribution (Q) is a very 
slow process. It can certainly not account for most of the numerical simulations 
and experiments of 2D turbulence and point vortex dynamics in which an equi- 
librium state is established extremely rapidly. This implies that a more violent 
relaxation mechanism must be at work in the system (see Sec. This is a 
remark of crucial importance because it means that all the results established 
previously, including the Boltzmann distribution (ph, must be revised. 



6 Violent relaxation of 2D vortices and stellar systems 
6.1 The Euler and the Vlasov equations 

In the preceding sections, we have focused our attention to the point vortex 
model as an idealization of more realistic flows which necessarily involve a contin- 
uous vorticity distribution. This approximation is interesting in a first approach 
because it leads to a system of N particles in interaction (like electric charges 
or stars) for which the methods of statistical mechanics are directly applicable. 
In addition, this model keeps the specificity of two-dimensional vorticity fiows 
such as long-range interactions between vortices and structure formation. How- 
ever, there are many different ways to approximate a continuous vorticity field 
by a cloud of point vortices and different approximations can lead to different 



statistical equilibrium states (this difficulty was underlined by Onsager |101|). 
Therefore, if we want to apply the results of statistical mechanics to realistic 
situations (e.g., geophysical flows) it is necessary to go beyond the point vortex 
model and develop a statistical mechanics for continuous vorticity fields. 

For flows of geophysical or astrophysical interest, the Reynolds numbers are 
so high that the molecular viscosity is not expected to play a crucial role in 
the dynamics. Therefore, these flows are described in the simplest model by the 
Euler-Poisson system 

^ + uVc^ = 0, (274) 
ot 

uj = -A^. (275) 

It can be recalled that these equations also model the early dynamics of a cloud 
of point vortices before correlations between vortices have developed (in that 
case, w is proportional to the one-body distribution function P(r,t), see Eq. 



(E33h). As discussed in Sec. 5.9, this is the regime of physical interest since the 



"collisional" relaxation of point vortices is in general very slow. 
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Similarly, for a majority of stellar systems, including the important class of 
elliptical galaxies, the relaxation time by two-body encounters is ~ 10^^ times 
larger than the age of the universe. Therefore, the dynamics of stars in a galaxy 
is essentially collisionless and appropriately described by the Vlasov-Poisson 
system 

% -4 <-) 

A$ = AttG j fd^v. (277) 

The morphological similarity of the Euler-Poisson and Vlasov-Poisson sys- 
tems is another manifestation of the close analogy between 2D vortices and stellar 
systems. If we make the correspondance between the vorticity and the distribu- 
tion function (uj <-!■ /) and between the stream function and the gravitational 
potential {i}) <P), these two equations describe the advection of a density by 
an incompressible flow with which it interacts via a Poisson equation. Then, the 
density is not advected passively by the flow but is coupled to its motion. This 
coupling is responsible for violent fluctuations of the stream function or gravita- 
tional potential. These fluctuations will mix the vorticity or the phase elements 
at small scales and induce a self-organization and the appearance of structures 
at larger scales (see Fig. This violently changing potential provides a mech- 
anism analogous to a relaxation in a gas, but the speciflcity of this relaxation 
is that it is collisionless and due to the long-range nature of the interactions. 
It is now clear that this "chaotic mixing" is the driving source of relaxation in 
two-dimensional turbulence and stellar dynamics. The kinetic theory presented 
in Sec. || for point vortices (and the one developed by Chandrasekhar for stars) 
is only valid in situations in which this chaotic mixing is prevented or has died 
away. 



6.2 The statistical equilibrium 

During the mixing process, the Vlasov-Poisson and the Euler-Poisson systems 
generate intermingled filaments at smaller and smaller scales. Therefore, a de- 
terministic description of the flow would require a rapidly increasing amount of 
information as time goes on. For that reason, it is appropriate to undertake a 
probabilistic description in order to smooth out the small scales and concentrate 
on the locally averaged quantities. This statistical approach, called the theory of 
"violent relaxation" , was introduced by Lynden-Bell |Q in 1967 for collisionless 
stellar systems and rediscovered independantly by Kuzmin [H, Miller ||95|,p6| 



and Robert & Sommeria [111| for two-dimensional vorticity flows. The analogy 
between these two statistical mechanics (including the relaxation towards equi- 
librium) was discussed in detail by Chavanis [^Qp^js^ , ^ . On the other hand, 
a rigourous justiflcation of this statistical approach has been given by Robert 
[ [109| by using the concept of Young measures and large deviations. This theory 
improves upon previous works based on a spectral representation of the flow 
[pO|, which do not respect all the conservation laws of the invisicd dynamics. In 
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Fig. 13. Violent relaxation of two-dimensional vorticity flows and self-gravitating 
systems. The left panel corresponds to the nonlinear developement of the Kelvin- 
Helmholtz instability. The equilibrium state is a large-scale vortex which is well- 
described by statistical mechanics [119|. The right panel corresponds to a simulation 
of the Vlasov-Poisson system showing a mixing process and the formation of a co- 
herent structure in phase space with a Fermi-Dirac distribution Q. This simulation 
is restricted to a one-dimensional system but the process remains the same in higher 
dimensions. 



the following, we present the statistical mechanics of violent relaxation for two- 
dimensional vorticity flows. A discussion of the statistical mechanics of violent 
relaxation in stellar systems, closely following the presentation of this paper, can 
be found in [||. 

In the statistical approach, the exact knowledge of the "fine-grained" or mi- 
croscopic vorticity field is replaced by the probability density p(r, a) of finding 
the vorticity level cr in r. The normalization condition yields at each point 

Jpir,a)da^l, (278) 

and the locally averaged (coarse-grained) vorticity is expressed in terms of the 
density probability in the form 

^= J pir,cT)ada. (279) 

More generally, the moments of the vorticity are defined by 

^ = J pir,cT)a''da. (280) 
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During the evolution, the energy 



E^l I ujiPcfr, (281) 



2 

is conserved as well as the total area of each level of vorticity 



7(a) = J p{r,a)d'v. (282) 

These last constraints are equivalent to the conservation of the Casimir integrals 
Ch = J h{uo)(fiY for any continuous function h. Such integrals are conserved by 
the Euler equation because the fluid particles keep their vorticity (on account of 
the transport equation cLu/dt — 0) and their surface (on account of the incom- 
pressibility of the flow). 

After a complex evolution the system is expected to be in the most probable, 
i.e. most mixed state, consistent with all the constraints imposed by the dynam- 
ics. We define the mixing entropy as the logarithm of the number of microscopic 
configurations associated with the same macroscopic state (characterized by the 
probability density p{r, a)). We divide the macrocells {r,r + dr) into h' microcells 
and denote by riij the number of microcells occupied by the vorticity level aj in 
the i-th macrocell. A simple combinatorial argument indicates that the number 
of microstates associated with the macrostate {uij} is 



W{{n,,})^l[N,\l[-^, (283) 

tl'ij ■ 



where Nj — is the total number of microcells occupied by cr, . W e have 

to add the normalization condition = equivalent to Eq. (278), which 



prevents overlapping of different vorticity levels. This constraint plays a role 
similar to the Pauli exclusion principle in quantum mechanics. Morphologically, 



the statistics (283) corresponds to a 4'^ type of statistics since the particles 
are distinguishable but subject to an exclusion principle ||90| ]. There is no such 
exclusion for point vortices (in the coUisional regime) since they are free a priori 
to approach each other without limitation. 

Taking the logarithm of W and passing to the continuum limit with the aid 
of the Stirling formula, we get 

S = - j cr) In p(r, a)d^vd(j. (284) 

The most probable macroscopic state is obtained by maximizing the mixing en- 



tropy (p84|) with fixed energy (281), global vorticity distribution (282) and local 



normalization (278). This problem is treated by introducing Lagrange multipli- 



ers, so that the first variations satisfy 



5S - (35E - J a{a)5-f{a)da - j Cir)5(^J p{r, a)da^ d^r = 0. (285) 
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The resulting optimal probability is a Gibbs state which can be expressed as 



where Ziijj) = e^(''^+^ and g{a) = e The normalization condition ( pTi 

leads to a value of the partition function Z of the form 



{a)e'^''^da, (287) 



and the locally averaged vorticity ( |279| ) is expressed as a function of ■0 according 
to 



This expression can be rewritten 

1 d\nZ 



OJ = — — - 



13 d-ip 



(289) 



Differentiating Eq. (288) with respect to ip^ we check that the variance of the 



vorticity can be written 



C^2=C^2_^2 (290) 



or, alternatively, 

1 d'^luZ , , 

Therefore, the slope of the function ZU = /(t/") is directly related to the variance 
of the vorticity distribution. Since loi > 0, we find that the function ZJ = f{ijj) 
is monotonic; it is decreasing for /? > and increasing for /3 < (it is constant 



for /3 = 0). Another proof of this result is given in |111 



Two particular cases are worth considering. If the local distribution of vor- 



ticity is Gaussian, then the uj — relationship (|288| ) is linear [£5 96|. On the 
other hand, in the case of a single level of vorticity (Tq (in addition to the level 
a = 0), the coarse-grained vorticity to = p(r, iJo)co takes explicitly the form 

uJ= ^V^- (292) 

This is formally similar to the Fermi-Dirac statistics. Here, the exclusion prin- 
ciple ZU < (Tq is due to the Liouville theorem (i.e., the conservation of the fine- 
grained vorticity) not to quantum mechanics. Because of the averaging pro- 
cedure, the coarse-grained vorticity can only decrease by internal mixing, as 
irrotational flow is incorporated into the patch do, and this results in an "effec- 
tive" exclusion principle. In the dilute limit ZJ ctq, Eq. (|292|) reduces to the 



Boltzmann distribution uj = Ae ^'^"'^ as in the point vortex model [|70|, 107 
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Similar results are obtained in the context of coUisionlcss stellar system de- 
scribed by the Vlasov-Poisson system . The equivalent of Eq. ( ^921 ) is the 
Lynden-Bell distribution function 

7= JVt^' (293) 

which formally coincides with the distribution function of the self-gravitating 
Fermi gas. In the non-degenerate limit / <ti ?7o, it reduces to the Maxwell- 
Boltzmann distribution / = Ae~^^°^^'^'^\ Therefore, the theory of violent re- 
laxation naturally explains the observed isothermal cores of elliptical galaxies 
without recourse to collisions M,B3| . 



6.3 The Maximum Entropy Production Principle 

Basically, a two-dimensional incompressible turbulent flow at high Reynolds 
numbers is described by the Euler-Poisson system. In principle, these equations 
determine completely the evolution of the vorticity field Lu(r,t). However, in 
practice, we are not interested by the finely striated structure of the fiow but 
only by its smoothed-out structure. Indeed, the observations and the numerical 
simulations are always realized with a finite resolution. Moreover, the coarse- 
grained vorticity ZJ is expected to converge towards an equilibrium state, the 



Gibbs state ( 286 ) , contrary to the exact vorticity field lo which develops smaller 
and smaller scales. There is also a technical difficulty to simulate an inviscid 
dynamics due precisely to the developement of this small-scale motion. Contour 
dynamics methods need to introduce a "surgery" and spectral codes a "viscos- 
ity" in order to prevent numerical instabilities. However, this artificial viscosity 
breaks the conservation laws of the Euler equations. What we would like to ob- 
tain is a set of relaxation equations which smooth out the small scales while 
conserving all the constraints of the Euler equation. Such a parametrization has 



been proposed by Robert & Sommeria [112| in terms of a phenomenological 
Maximum Entropy Production Principle (MEPP). 

Let us decompose the vorticity uj and velocity u into a mean and a fiuctuating 
part, namely u = oJ+uj, u = u-l-u. Taking the local average of the Euler equation 



( P74D , we get 

|^ + V(cJu) = -VJ^, (294) 

where the curr ent = uiu represents the correlations of the fine-grained fluc- 
tuations. Eq. ( |294| ) can be viewed as a local conservation law for the circula- 
tion r = J uj(Pr. To apply the MEPP, we need to consider not only the lo- 
cally averaged vorticity field U but the whole probability distribution p(r, cr, t) 
now evolving with time t. The conservation of the global vorticity distribution 
7(ct) = / p(Pr can be written in the local form 

|^+V(pu) = -VJ, (295) 
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where J(r, a, t) is the (unknown) current associated with the vorticity level a. In- 
tegrating Eq. ( |295D over all the vorticity levels cr, using Eq. (278), and comparing 
with Eq. (P), we find the constraint 



J J{r, a,t)da ^ 0. 



(296) 



Multiplying Eq. ( 295 ) by a, integrating over all the vorticity levels, using Eq. 
( ^79|) , and comparing with Eq. (294), we get / 3{r,a,t)ada = Jui- 



We can express the time variation of energy in terms of J, using Eqs. (281) 
and (294), leading to the energy conservation constraint 



E = / J^VV- d^r = 0. 



(297) 



Using Eqs. (^84|) and (295), we similarly express the rate of entropy production 
as 

3W{lnp)d'^rda. (298) 



5* 



The Maximum Entropy Production Principle (MEPP) consists in choosing 
the current J which maximizes the rate of entropy production S respecting 
the constraints E = 0, J Jda = and J ^da < C{r,t). The last constraint 
expresses a bound (unknown) on the value of the diffusion current. Convexity 
arguments justify that this bound is always reached so that the inequality can 
be replaced by an equality. The corresponding condition on first variations can 
be written at each time t: 



6S - (3{t)6E - 



and leads to a current of the form 



J = -D{r,t) 



(299) 



(300) 



The Lagrange multiplier ^ has been eliminated, using the condition (|296|) of 



local normalization conservation. The conservation of energy (297) at any time 



determines the evolution of the Lagrange multiplier (3{t) according to 



J DViJV-ipd'^r 



The entropy production (|298| ) can be written 



(301) 



-(Vp + Ppia - uj)\7il))d^Yda + (3 I Jia - UW^d^rda. (302) 
P 
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Using Eqs. (296) and (297), the second integral is seen to cancel out. Inserting 
Eq. (|300[) in the first integral, we find 



S = / —(fvda, 
Dp 



(303) 



which is positive provided that I? > 0. A stationary solution S' = is such that 
J = yielding, together with Eq. (|300|), 



For any reference vorticity level ctq, it writes 



0. 



(304) 



(305) 



Substracting Eqs. (|304| ) and ( |305| ) , we obtain Vln((o//9o) + I3{(J 
which is immediately integrated into 



cr) 



1 



Z{v) 



ao)V^ = 0, 



(306) 



where ^-^(r) = po(r)e'^'"°'^('") and g{a) 



, A{a) being a constant of inte- 



gration. Therefore, entropy increases until the Gibbs state (^86|) is reached, with 
P = limj^oo f3{t). Furthermore, we can show that a stationary solution of these 
relaxation equations is linearly stable if, and only if, it is an entropy maximum 
(in preparation). Therefore, this numerical algorithm selects the maxima (and 
not the minima or the saddle points) among all critical points of entropy. When 
several entropy maxima subsist for the same values of the constraints, the choice 
of equilibrium depends on a complicated notion of "basin of attraction" and not 
simply whether the solution is a local or a global entropy maximum (see Ref. 
[p6| in a related context). 



The relaxation equations (295), ( [300 ) and (301) can be simplified in the single 
level approximation. In that case, Eq. (294) is explicitly given by 



duj 
'dt 



(307) 



In the dilute limit oJ <C ctq, it takes a form similar to the Fokker-Planck equation 
( p06D obtained for point vortices. These equations both involve a diffusion and 
a drift, but these terms have a different physical interpretation in each case. 

It is also instructive to apply this thermodynamical approach to the Vlasov- 
Poisson system. In the single level approximationjjt leads to the following equa- 
tion for the coarse-grained distribution function ||5^ , |3T| : 



df df Bf d 
— + V— + F— = — 
dt dr dv dv 



D\^ + m.f{m-f)^ 



(308) 



which is a generalized form of the familiar Fokker-Planck equation (|20^ ) recov- 
ered in the non degenerate limit / <C rj^. W e can check that Eq. (308) returns 
the Lyndcn-BcU distribution function (293) at equilibrium. 
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The diffusion coefficient D is not determined by tlie MEPP as it depends on 
the unknown bound C on the current. For the purpose of reaching the Gibbs 
state ( ^86[ ) , the diffusion coefficient can simply be chosen arbitrarily (but with a 
positive value in order to ensure entropy increase). However, the precise form of 
the diffusion coefficient is important in order to determine the relaxation time 
and take into account kinetic confinement and incomplete relaxation (see below). 
In the context of 2D turbulence, Robert & Rosier [110| have proposed a simple 
evaluation of D by using an analogy with the diffusion of a passive scalar uj 
subjected to a turbulent velocity field u = u-|- u (see also |54 , Appendix B). In 
that case, the mean value tU satisfies a convection-diffusion equation 



'dt 



(uV)cJ= V(L»VtU), 



with a diffusion coefficient given by 



(309) 



(310) 



where r is the decorrelation time of the velocity fluctuations. Equation (307) 
reduces to Eq. (31C) when /3 = 0, i.e. when the energy constraint is not active. 
In the Euler-Poisson system, the velocity is produced by the vorticity itself via 
the Biot & Savart formula 



u(r,t)= J Lj{r',t)V{r' ^ r)(fr', 



(311) 



where V(r' ^ r) = -(l/27r)z x (r' - r)/|r' - rp. Inserting Eq. (|31l| ) in Eq. 
(310), we obtain 



D = ^J Cj{r',t)u;{r",t)V{r' r)V(r" -> r)d^r'd^r" 



(312) 



If we neglect the spatial correlations of the vorticity fluctuations on scales larger 
than e, the resolution scale, we have 



Co{r',t)Cj{r", t) = eV(r',t),5(r' - r"), 



and Eq. (312) reduces to 



D 



\ c;)2(r',t)V^(r' ^r)d^r'. 



(313) 



(314) 



Making a local approximation and introducing an upper cut-off a, we obtain 

D = -/Cj^iv, t) r f^] 2<de, (315) 



with J)2 
coefficient 



UJ . This leads to the following expression for the diffusion 



D = (cj2 - uj' 



-2^^1n 



(316) 
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We note that the diffusion coefficient vanishes in regions where there is no fluc- 
tuation of the vorticity at small scales, i.e. ljJ2 = 0. At the contact with the 
unmixed flow, the diffusion current also vanishes resulting in a conflnement of 
the vorticity. This leads to the concept of incomplete relaxation: in a large-scale 
turbulent flow, the Gibbs state ( |286| ) is satisfled only in restricted regions of 
space where mixing is sufficiently efficient to justify an ergodic hypothesis. Out- 
side these domains, the relaxation is slowed down or even stopped. This kinetic 



confinement is illustrated by the numerical simulations of Robert & Rosier | 110 | 
and further discussed by Chavanis & Sommeria [si}] . In this viewpoint, the vor- 
tices of two-dimensional turbulence can be considered as restricted equilibrium 
states or maximum entropy bubbles separated from each other by an almost 
irrotational background. A similar kinetic confinement can be advocated in the 
case of stellar systems |^,^ to account for incomplete relaxation Q and solve 
the infinite mass problem. 

6.4 Recent developements 

An interesting problem in fluid dynamics is to obtain a classiflcation of the 
"zoology" of vortices (monopoles, translating and rotating dipoles, tripoles...) 
met in two-dimensional flows. The statistical mechanics approach presented in 



Sec. 6.2 provides a general framework to tackle this problem as it selects the 
most probable structures among all possible solutions of the 2D Euler equations. 
However, the prediction is not straightforward because, in the general case, we 
have to take into account an inflnite set of constraints, namely the conservation of 
all the Casimirs C'h, in addition to energy E, angular momentum L and impulse 



P. A numerical algorithm has been developed by Turkington & Whitaker [124| to 
solve this problem and several calculations have been performed in rectangular 
or circular domains for a flnite number of vorticity levels and for particular values 
of the integral constraints. However, many structures are found and it is difficult 
to have a clear picture of the bifurcation diagram in parameter space. For that 
reason, Chavanis & Sommeria |49y5l[] have considered a particular limit of the 
statistical theory, the so-called "strong mixing limit" , in which the study of these 
bifurcations can be performed analytically. This limit corresponds to (3a^ <^ 1 
so that the equations of the problem can be expanded in terms of this small 
parameter (this is like the Debye-Hiickel approximation in plasma physics). To 
zeroth order in the expansion, the density probability of each level is uniform 
which corresponds to a completely mixed state. To first order, the relationship 
between vorticity and streamfunction is linear and this can justify an inviscid 
minimum enstrophy principle (for the coarse-grained enstrophy ^' — JuJ^cPr) 
In that case, the equilibrium flow only depends on E, L, P and the flrst 
moment P of the vorticity (the flne-grained enstrophy P^'^' = J cu'^cPr serves 
as a normalization factor) . This particular limit of the statistical theory already 
exhibits a rich bifurcation diagram and often provides a good approximation 
of more general situations (in particular for weakly energetic flows). It is also 
possible, in principle, to go to higher orders in the expansion in which case more 
and more vorticity moments 13,14,... are necessary to describe the structure 
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of equilibrium. This limit makes therefore a hierarchy between the constraints 
as it shows that, in many situations, only the lowest moments of the vorticity 
are important to characterize the equilibrium state. Because of this hierarchy, 
we can make some relevant predictions without the complete knowledge of the 
initial condition. 





Fig. 14. Jupiter's great red spot 



The statistical mechanics approach can be extended immediately to the 
quasi-geostrophic (QG) equations |94 by simply replacing the vorticity by the 
potential vorticity (PV). The formalism has also been generalized by Chava- 
nis & Sommeria to the shallow-water (SW) equations. These equat ions are 
more relevant to describe geophysical flows than the 2D Euler equations 1104]. In 
the QG approximation, Bouchet & Sommeria fl^ have explained the formation 
of jets and vort ices in planetary atmospheres in terms of statistical mechanics, 
as initiated in [ 118 |. In particular, the annular jet structure of Jupiter's Great 
Red Spot (see Figs. |l4-15) is reproduced and explained as the coexistence of 
two thermodynamical phases in contact (a picture which is rigorously valid in 
a small Rossby radius expansion). These results can be extended to the more 
general case of shallow water equations ||l^. In this geophysical context, the de- 
formation of the fluid surface tends to reduce the range of interaction between 
vortices. This situation is comparable to what happens in a neutral plasma: the 
interaction between vortices is shielded on a distance of the order of the Rossby 
length, the analogous of the Debye length in plasma physics. 

The relaxation equations presented in Sec. 3.2 provide a convenient param- 
eterization of sub-grid scale eddies which drives the system toward statistical 
equilibrium by a continuous time evolution. Such relaxation equations can be 
used both as a realistic coarse resolution model of the turbulent evolution, and 
as a method of determination of the statistical equilibrium resulting from this 
evolution. However, these equations do not preserve the invariance properties of 
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Fig. 15. Relaxation towards statistical equilibrium in a QG model of Jupiter's great red 
spot (from Ref. jl^]). Three successive potential vorticity (PV) fields are represented 
as grey levels. The initial condition (top) is made of small PV patches. These patches 
organize into vortices (middle) that eventually merge into a single one (bo ttom) . This 
sequence is obtained with the relaxation equations described in Sec. 6.3, applied to 
the Q.G. situation: entropy increases with time while energy is exactly conserved. 
At equilibrium, the vortex is an oval spot of quasiuniform PV surrounded by strong 
gradients, corresponding to an annular jet. 



the Euler equations by a change of reference frame. In addition, the conservation 
of energy is enforced by a formal Lagrange multipher P{t) which is uniform in 
space. This may be a limitation to describe large-scale turbulent flows which 
organize locally in several types of structures with a different temperature. A 
generalization of the MEPP has been proposed by Chavanis & Sommeria 
as an attempt to solve these difficulties. More general, but also more complex, 
relaxation equations arc obtained. They involve a space dependant temperature 
which tends to be uniform in each vortex but with an a priori different value 
from one vortex to the other. A simplified version of these relaxation equations 
has been solved numerically by Kazantzev et al. [ [78| in an oceanic context. 

Despite its practical interest, the main drawback of the MEPP is it s ad 
hoc nature. In we have attempted to justify the relaxation equation (307) 
from first principles starting directly from the 2D Euler equation. A systematic 
derivation can be obtained in the so-called quasilinear approximation. This ap- 
proximation is well-known in plasma physics and stellar dynamics ]71,114,42[ 
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for the Vlasov-Poisson system and we tried to extend it to the Euler-Poisson 
system. Substracting Eq. (p7j 



and (294) and neglecting the non hnear terms in 
the equation for cj, the basic equations of the quasihnear theory are 



'dt 



uVuj = — VcDu, 



— + uVw 

at 



-uVw. 



(317) 



(318) 



It is possible to solve Eq. ( pl8D formally with the aid of G reen 's functions and 
substitute the resulting expression for Cb{v, t) back into Eq. ( |317| ). Implementing 
a closure relation of the form (313), a kinetic equation can be obtained for tU. In 
the single level approximation, it reads p%: 



at 



ds / (^Y'V^'{Y' ^ r)t 



^r^^ Jq 

X ^ V^{r' ^ rP'ia, - -')^ + V^r ^ r')^(ao - -)^ 



(319) 



Assuming in addition that the decorrelation time t is short, as in the stochas- 



tic model of Robert & Rosier |11C], we are led to the following equation 



d 



K^^ir' - r) = 



(320) 
(321) 



with ^ = r' — r. This equation includes a diffusion and a drift, as in the MEPP, 
but these terms are obtained here directly from a local average of the Euler 
equation. On the other hand, the conservation of angular momentum results from 
the symmetry of the diffusion current instead of an ad hoc Lagrange multiplier. 
This symmetric structure respects in addition the invariance properties of the 
Euler equation. Finally, a if-theorem for the Fermi-Dirac entropy 



1 In 1 



(322) 



can be derived from this kinetic equation instead of being postulated as in the 
MEPP. Our approach provides therefore an alternative, dynamical, justification 
of the mixing entropy introduced by Miller ||9^ and Robert & Sommcria [112| 
at statistical equilibrium 



Equations (B20) and (321) are similar to the kinetic equations (23S) and ( p3£ 
obtained for point vortices. There are, however, two important differences: (i) 
the drift and the diffusion involve the product ZU x (oq — u) instead of (w) . This 
nonlinearity ensures that the constraint ZU < ctq is satisfied at any time, (ii) The 
diffusion coefficient is proportional to the circulation croe^ of a completely filled 



2D vortices and stellar systems 



71 



macrocell, instead of the circulation 7 of a point vortex. In general, uoe^ 3> 7 so 
that the relaxation by coUisionless mixing is much more rapid than the coUisional 
relaxation. From the above theory, we find that the time scale of the violent 
relaxation is of order to, the dynamical time, whereas the collisional relaxation 
of point vortices operates on a time scale ^ {N/ In A^)^^). 

If we are close to equilibrium, we can implement a "thermal bath approxi- 
mation" an d rep lace the vorticity ui' = uj{ r' , t) by its equilibriimr form ( ^92| ) [3^ . 
Then, Eq. (|20|) reduces to the equation ( ^ ) derived from the MEPP and the 
diffusion coefficient coincides with the estimate ( ^16[ ) based on the passive scalar 
model (cj2 — ^i^o ~ ^) in the single level approximation) . All these results are 
satisfactory. However, Eq. (32C) does not conserve energy and this marks a flaw 
in the previous description. It is probable that this constraint is broken by the 
small correlation time hypothesis, but we do not see how to simplify Eq. ( |319| ) 
further without this assumption. We could try to satisfy the energy constraint 
by introducing a (5-function term in Eq. ( |320| ) as in Eq. ( |245| ) but it is not clear 
how one can justify this procedure in the present context. The conservation of 
energy (which is intimately related to the form of the drift term) is the most 
serious problem that we have encountered in trying to develop a kinetic theory 
of 2D turbulence. Up to date, it has not been answered satisfactory. Further 
scrutinity should be given to the non markovian kinetic equation (319) which 
may conserve energy although we were not able to prove it. 

It has to be emphasized that the quasilinear theory is not a theory of violent 
relaxation in the usual sense since it only applies to the late quiescent stages 
of the relaxation ("gentle relaxation"), when the fluctuations have weaken and 
a linearization procedure can be implemented. In order to give a more relevant 
description of violent relaxation, it will be important in future works to make 
a link between statistical mechanics and chaotic dynamics, in particular in the 
point vortex model. This may give a new estimate of the diffusion coefficient and 
of the decorrelation time r, which in this context could be related to a Lyapunov 
exponent. We feel that this track is an important one to make progress in the 
understanding of 2D turbulence and point vortex dynamics. 

In the quasilinear theory, it is implicitly assumed that the decomposition 
uj = uj + Lu is obvious and that uj should be regarded as a statistical average 
of CL'(r,t,C) over different realizations C of the flow. This implies in particular 
that ZJ — UJ. However, a different approach is considered by Laval et al. and 
Bouchet ll^ who define To as the convolution of cj(r, t) with a Gaussian window 
of size e (or b y a t runcation in Fourier space) . In that case uj ^lJ and new terms 
arise in Eq. ( 317 ), in particular a term — ufa — cJ u which dominates over 
the others [[l2[. When e — > 0, this term becomes equivalent to an anisotropic 
diffusion e'^S'VlU with a diffusion coefficient (or turbulent viscosity) related to 
the stress tensor ^7^'^- = diuj. This term alone conserves energy but there exist 
directions in which the viscosity is negative leading to instabilities. In order to 
circumvent this difficulty, Bouchet [|l2j proposes to project the diffusion current 
on directions in which the viscosity is positive and to introduce a drift term, as in 
the MEPP, in order to recover the conservation of energy lost by this procedure. 
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This leads to an operational subgridscale model of 2D turbulence which appears 
to be more efficient than other parametrizations. As for the MEPP, its drawback 
is its ad hoc nature but it is not clear at present if it will be possible one day to do 
much better and derive a parametrization of 2D turbulence from first principles 
as attempted in the quasilinear theory. 



6.5 The limits t — > +00 and N +00 

We have indicated previously that a system of point vortices or point mass stars 
can achieve two successive equilibrium states. On a short time scale, the corre- 
lations between particles have not yet developed and the system is described by 
the Vlasov (or Euler) equation. In this regime, the dynamics is collisionless. Yet, 
for systems with long-range interactions, the collective nature of the evolution 
is responsible for an effective relaxation process, called violent relaxation, which 



leads to a metaequilibrium state (292) (293) on a very short time scale. On a 



longer time scale, the fluctuations of the potential have died away and the devel- 
opement of correlations between stars or between point vortices leads to another, 
slower, relaxation process. This corresponds to the "coUisional" regime. This sec- 
ond process is more standard and leads to a true equilibrium state (^)( |lOl| ). 
In the case of continuous vorticity fields (instead of point vortices) , the second 
stage is replaced by a viscous decay of the vorticity due to inherent viscosity. 

We can discuss these two successive equilibrium states in a slightly different 
manner. Let us consider a fixed interval of time and let the number of particles 
N +00. In that case, the system is rigorously described by the Vlasov (or 
Euler) equation, and a metaequilibrium state is achieved on a timescale inde- 
pendant on N. Alternatively, if we fix N and let t +00, the system will relax 
to the true equilibrium state resulting from a coUisional evolution. These two 
equilibrium states are of course physically distinct. This implies that the order 
of the limits N +00 and t — > +00 is not interchangeable. 

In the process of violent relaxation, the statistical mechanics is not as firmly 
established as in the coUisional regime, although it is often the process of most 
interest. Indeed, the mixing required for the validity of the ergodic hypothesis 
is fed by the fluctuations of the potential. As these fluctuations decay as we ap- 
proach equilibrium (by deflnition!), the mixing becomes less and less efficient and 
this can lead to an incomplete relaxation. Since the Boltzmann-Gibbs entropy 
does not always give a good description of the equilibrium state, it has been 
proposed sometimes to use a wider class of functionals to describe the process 



of violent relaxation |122|. Among them, the g-entropies introduced by Tsallis 
[ |123| have been shown to give in some cases a good fit of the equilibrium state. 
However, since the value of q needs to be adjusted in each case, it is not clear 
whether this agreement is the signal of a generalized thermodynamics or just 
a coincidence n. It is clear that the Boltzmann entropy does not always give a 



^ In 2D turbulence, Boghosian justifies a form of minimum enstrophy principle 
from Tsallis thermodynamics in order to interpret the experimenal results of Huang 
& Driscoll in a magnetized plasma. This is because the enstrophy /2 ~ J Lj^d^r 
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perfect description of the equilibrium state but there is no convincing reason, 
up to date, why the system would select another "universal" form of entropy in 
the context of violent relaxation Q|. In any case, the formalism developed by 
Tsallis and co-workers is nice to generalize to a wider class of functionals the 
results obtained with the Boltzmann entropy. Since this generalization is often 
analytically tractable (as it leads to power laws), this may explain the interest 
and the attractive nature of this approach. 



7 Conclusion 

The statistical mechanics of 2D vortices and stellar systems appear to be re- 
markably similar despite the different nature of these systems. We have tried 
to develop this analogy in different directions. First of all, the iV-star and A^- 
vortex problems both involve an unshielded long-range potential generated by 
the density of particles themselves. At statistical equilibrium, these systems are 
described by a Boltzmann-Poisson equation whose solutions characterize orga- 
nized states (at negative temperatures for vortices). In the case of stars, the 
relaxation towards equilibrium can be viewed as a competition between a diffu- 
sion and a friction. We have proposed to describe the relaxation of point vortices 
similarly in terms of a diffusion and a drift. The diffusion is due to the fluctua- 
tions of the force experienced by a star or to the fluctuations of the velocity field 
moving a vortex. The statistics of these fluctuations can be studied by similar 
mathematical methods. The fluctuations of the gravitational field are described 
by a particular Levy law, called the Holtzmark distribution, and the fluctuations 
of the velocity of vortices are described by a marginal Gaussian distribution, in- 
termediate between Gaussian and Levy laws. The friction experienced by a star 
is due fundamentally to the inhomogeneity of the velocity distribution of the 
stellar cloud. Analogously, the drift experienced by a vortex results from the 
spatial inhomogeneity of the vortex cloud. The friction and the drift can be un- 
derstood similarly in terms of a polarization process and a back reaction of the 
system. In the thermal bath approximation, the coefficients of friction and drift 
are given by an Einstein relation and the one-body distribution function sat- 
isfies a Fokker-Planck equation. Further away from equilibrium, the collisional 
dynamics of stars is described by the gravitational Landau equation. We have 
derived a new kinetic equation that should be appropriate to the "collisional" 
relaxation of point vortices. 

A system of stars or vortices can also undergo a form of violent relaxation. 
This is essentially a collisionless process driven by the rapid fiuctuations of the 
potential as a result of collective effects (chaotic mixing). In this collisionless 



is a particular q-entropy. However, this is essentially coincidental and the minimum 
enstrophy principle can lead to inconsistenc ies a s discussed in lid]. In astrophysics. 



Tsallis entropies lead to pure polytropes [106]. These distribution functions are 
known for a long time but they do not give a particularly good description of ellip- 
tical galaxies or other stellar systems. Therefore, the relevance of Tsallis generalized 
thermodynamics in 2D turbulence and stellar dynamics remains questionable lid 
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regime, the stars and the vortices are described by the Vlasov-Poisson and the 
Euler-Poisson systems. These equations present a similar structure and a statis- 
tical mechanics can be developed to predict the "most probable state" resulting 
from a complex evolution driven by a mixing process. The relaxation towards 
equilibrium can be incomplete and an out-of-equilibrium study is necessary to 
understand what limits relaxation and causes a kinetic confinement of the sys- 
tem in a "maximum entropy bubble" . Relaxation equations have been derived 
either from a heuristic Maximum Entropy Production Principle or from a more 
controllable kinetic theory, in an asymptotic regime of the dynamics (gentle 
relaxation) in which a quasilinear approximation can be implemented. 

It should be noted that many results presented in this paper, in particular 
those corresponding to the kinetic theory of 2D vortices developed in Sees. ||-||, 
are very recent and need to be completed and further discussed. In particular, it 
appears indispensable to carry out extensive numerical simulations to test their 
relevance and determine their domains of applicability. It is plausible that the 
true dynamics of stars and vortices is more complex than the picture that has 
been given here. In addition, the description of chaos in these systems has not 
been addressed at all in this paper although it is presumably an essential ingre- 
dient to understand their dynamics. We are thus far from reaching a complete 
understanding of these systems with long-range interactions. We feel, however, 
that the analogy between the statistical mechanics of stars and vortices that we 
have investigated is correct in its mains lines and should lead again to fruitful 
developements. 



Appendix A: The calculation of the diffusion coefficient 

In this Appendix, we calculate the diffusion coefficient of a point vortex evolving 
in an inhomogeneous vortex cloud, using the Kubo formula. To evaluate the 
velocity autocorrelation function C(t) — {V{t)V{t — r)), we shall assume that 
between t and t — t, the point vortices follow the streamlines of the equilibrium 
flow. This is a reasonable approximation in the case of strong shears. 

7.1 Unidirectional flow 

The trajectory of a point vortex advected by a unidirectional equilibrium flow 
is simply: 

y{t-T)=y{t), (323) 
x{t~T)^x{t)~{V)MT. (324) 



The velocity auto-correlation function appearing in Eq. (185) can be written 
explicitly 

^(^) / dxidyi- T^— T ^2 W 

47r^ J (xi - xy + (yi - yy 

Xi — X 



{xi - xy + {yi - yY 



{t-T)P,g{y), (325) 
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where we have used Eq. (p^. The second term involves the quantity 

(xi - x){t _ r) = - x + {{V),M - {V),q{y)y- (326) 

Since the integral in Eq. ( |325| ) diverges as ri — > r, we can make a local approx- 
imation and expand the velocity difference in a Taylor series m. yi — y. To first 
order, we have 

{V)eM - {V)eM ^ -S{y){y, - y), (327) 

where S{y) is the local shear of the flow. Introducing the variables X = xi — x 
and Y = yi — y, we obtain 

The integration over X can be performed easily since the integrand is just a 
rational function of polynomials. After straightforward calculations, we find 

The integral over Y diverges logarithmically for both small and large Y. The 
reason for this divergence has been explained in Sec. |^. Introducing two cut-offs 
at scales d (the inter-vortex distance) and R (the domain size) and noting that 
\n{R/d) ~ i In A^, we obtain 

Cir) ^^\nN—^l-—PM- (330) 
For r +00, the correlation function decreases like r^^. This is a slow decay 



but it is sufficient to ensure the convergence of the diffusion coefficient (185). 
Using 

/■* N-y'^ In TV /I \ 

C{r)dr = _^^arctan(^-|i;(y)|ijPe,(y), (331) 



and taking the limit t — > -)-oo, we establish Eq. ( 189 ) 



7.2 Axisymmetric flovir 

In an axisymmetric flow, the trajectory of a point vortex takes the simple form: 

r{t - r) = r(i), (332) 



e{t -t)= e{t) - MfiMr. (333) 
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As indicated in Sec. 4.2, we are particularly interested by the r{t)r{t — r) com- 
ponent of the correlation tensor. Let us introduce the separation (5r = ri — r 
between the field vortex f and the test vortex. In the local approximation, 5y 
can be considered as a small quantity. Therefore, we can write 



(334) 



(fvi = (P{5v) = dXdY. 



(335) 



With these notations, the correlation function appearing in Eq. ( |188D can be 
rewritten 



C(r) 



47r2 



Peg(r) j dXdY 



(336) 



Now, 



and 



Y{t -t) = Sr(t ~ t) = ri{t - t) - r{t - t) 
= ri{t)-r{t)^Y{t)^Y, 



X{t - r) = r{t - T)5dlt -t)= r{t - T){9i{t - r) - 9{t - r)) 
\ \ n r 



(337) 



(338) 



In the local approximation, we can expand the last term in Eq. ( |338| ) in a Taylor 
series in ri — r. This yields 



X{t^r)=r{e,{t)-e{t))-r±(^^^^ 

dr \ r 



(ri -r)T = X - S{r)YT, 



(339) 



wher e S (r) is the local shear of the flow. Substituting Eqs. (337) and ( |339|) in 
Eq. (pel), we get 



C(r) 



4^ 



Peg(r) J dXdY 



X 



X - E{r)YT 



X^ + Y^{X- S{r)YT)^ + y2 • 



(340) 



This integral is similar to Eq. (328), so we again obtain Eq. (n8 



Appendix B: Calculation of the memory function 



In this Appendix, we calculate the memory function that occurs in Eq. ( |236| ). If 
we assume that P = P{r,t), then Eq. (|236|) simplifies to 
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where 



Jr = -N 



dT J d^riVrit)il ^ 0)t 
( BP BP 1 



(342) 



and where V^(t)(l 0) is the coniponeiit of the vector V(l 0) in the direction 
of r(t). If we denote by {r{t),6{t)) and (ri(t), 0i(t)) the polar coordinates that 
specify the position of the point vortices and 1 at time t, we easily find that 



v;(i)(i^o) = -, 



ri sin( 



27r rf + 7-2 — 2rri cos( 



(343) 



We shall assume that between t and t — t, the point vortices follow circular 
trajectories with angular velocity f2{r, t). In that case, r{t — T) = r and 9{t — T) = 
6 — n{r, t)T. Then, we obtain 



v;(t-r)(i^o) 



7 



ri sin(6l -6^,- AQt) 



27r ri 



2rri cos{0 - 0i - AQt) 



with 



= f2{r,t) - Q{ri,t). 



(344) 



(345) 



We find similarly that Vr^(t-T){^ ^ 0) = ■^Vr(t-T){^ ^ 0)- Our previous as- 
sumptions also imply that P{r{t — T),t — t) ~ P{r,t) between t and t — t. In 
words, this means that the correlation time is smaller than the time scale on 
which the average vorticity changes appreciably. We do not assume that it is 
much smaller as in Sec. 5.5, so this approximation is not over restrictive. In that 



case, the diffusion current becomes 



Jr = -N 



+ 00 



dT 



2ir 



dOi 







+ 00 







rridriK(t)(l ~^ 0) 



xK(*-r)(l->0) 



1 BP 1 BPi 

-^17^ ^TT^ 

r or ri ori 



(346) 



where the time integral has been extended to +00. We now need to evaluate the 
memory function 



M 



+ 00 



dT 



2ir 



rf^^iv;(t)(i^o)K(*-r)(i->o). 



Introducing the notations (j) — 61 



- 6 and 
2rri 



<1, 



(347) 



(348) 
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we have explicitly 



M ^ { —] / dr d4> -— -4—. (349) 

K^irr) Jo Jo 1-ACOS0 1 - \cos{(j) + AQt) ^ ' 



This can also be written 



M 

where 



(i^j I '^^ j d(bV\(j))V'{(b + Af2T), (350) 



^(0) = ln(l - Acos0). (351) 
We now write the function V{(j)) in the form of a Fourier series, 



+ 00 „^ 

¥{<!>) = ^ a„e™^ with «n = 2^ / y(</>)e""*c 



(352) 



The memory function becomes 



A/ 



= --f^j y dr J del) J2 "ma«a„,e'("+'")"^e*™'^^^ (353) 



n.m— — oo 



Carrying out the integrations on and r using the integral representation of the 
delta function 



1 



+ 00 



Six) - — / e~^P-dp, (354) 



we are left with 



2 j-oo 2 +og 

^^"""g^ X! ?ima„am(5„ _„(5(mZir2) = — (5(Zil2) ^ |n|a^. (355) 

n,m— — oo n— — oo 

It remains for us to evaluate the series that appears in the last expression of the 
memory function. Using the identities 



ln(l - Acos^!))cos(n0)d0 = -^^i - - 1 ) (n > 0), (356) 



^ ln(l-Acos0)d(/) = 7rlnQ + ^-^^^^, (357) 



and the definition (34S) of A, we find that ao < oo and, for n > 0, 

(358) 



2rri J n\r 



> 
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where r> (resp. r<) is the biggest (resp. sniallest) of r and ri. Therefore, the 
value of the series is 



+ OC 



+ OC 



n — — oo n— 1 n— 1 ^ 

The memory function takes the form 

„2 



2n 



-2 In 



2n 



M 



7 

4j,2 



(5(zir2) In 





/ \ 2n 


1 - 


V 







and the diffusion current in the axisymmetric case can be written 



(359) 



(360) 



4r 



2 r+oo 



ridriS(f2 — i7i) In 



1-1 — 

r> 



r dr Ti dri 



(361) 



This leads to the kinetic equation (243). 
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